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1 Introduction 

The well known theorems of Hawking and Penrose Q proof the existence of spacetime singularities in a broad 
class of GR solutions under certain mild conditions on the matter energy-momentum tensor. Well known 
examples are black holes and the initial cosmological singularity in the standard Big Bang scenario. It is 
commonly believed that in this large curvature regions of the gravitational fields GR fails to be a completely 
consistent description and should be superseded by the more fundamental theory, in all probability, quantum 
in its nature. At present the most promising candidate for the consistent quantum gravity theory and for a 
unification of fundamental interactions are string theory and its recent extensions such as M-theory In 
these theories the characteristic scale in which the stringy effects become important is Planckean size which is 
much smaller than the scale we can probe in recent high energy experiments. Cosmology is only available way 
to test the consequences of the theory, since the cosmological predictions originate from physics of the early 
Universe, when space-time curvature may have been of Planckean strength (about the possibilities the string 
scale being close to the electroweak scale see [|)-[||). The today detectable remnants of primordial processes 
could be a test to which to compare the theory and a number of problems of standard cosmology, like singularity, 
dark matter and large scale structure, could be solved. 

The standard approach to string cosmology is based on the analysis of time dependent solutions of the 
effective field theory equations (see, for example, refs The effective theory of gravity resulting from 

string theory at low energies includes important modifications of GR due to the presence of the additional degrees 
of freedom such as dilaton, axion and etc. These fields couple to each other and to gravity nonminimally, and 
can influence the dynamics significantly. The study of cosmological consequences of these modifications has been 
an area of much active research in the past. Many authors have considered various aspects of the solutions, 
both in the E-frame and in the string frame. The simplest case is the pure gravi- dilaton solution Q, Q, (m), 
p7{ , p4[ , p5| , p7y , which is used in to construct the pre-big bang scenario of cosmological inflation. One 
can obtain further generalizations of this solution for the models with spatial curvature |17[-|27| and with the 



antisymmetric tensor fields of both NS-NS and R-R types |lj, [£7|, ||-[^,J^,|£g,|]^,Q,[|B5|,|7g] (and 
references therein), either by directly solving the corresponding cosmological equations or by using generating 
techniques. Such techniques include various dualities (see, e.g., pig , p5[ , [^0|-|6^|, j7lj ) and a dimensional 
reduction of the higher dimensional black hole solutions |^6| , |58| . Much of the emphasis in previous work on 
string cosmology has been on the problem of cosmological singularity. The attempts to address it included 
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the use of winding modes |13|, [ p2| , higher-derivative corrections |34f , |3q l-|4(|, |p0| , models with generalized 
scalar-tensor couplings |jS3f, |67|, inhomogeneous phase [pi] , |59fl , [f75f and etc. 

A large number of cosmological models have been derived in the past from a variety of vacua of M-theory 



|1[,|^,|^,J72[,I7§,179|. These include solutions of type II strings with R-R and NS-NS background fields. 
In particular, it has been pointed out that cosmological solutions can be generated from p-brane solutions by 
inverting the roles of the time and radial spatial coordinate. A new way of investigating string cosmology is 
connected to the conjecture due to Maldacena |8]{] , which relates string theory in Anti-de-Sitter spacetime to a 
conformal field theory living on its boundary (see ||S2||). 

One of the basic features of superstring theories is the existence of scalar fields, called moduli, which couple 
with gravitational strength and are massless to all orders in perturbation theory. The most general modulus field 
is the dilaton, whereas the other moduli parametrize the size and the shape of the extra compact dimensions. 
These fields are natural partners of the metric and play an important role in cosmology. A number of parameters 
in the string effective action, such as gauge and gravitational couplings are determined by expectation values 
of the moduli. As a consequence the cosmological variations of moduli fields will lead to the corresponding 
variations of the physical constants. To avoid conflict with observations conventionally it is assumed that non- 
perturbative phenomena create a potential which generate masses for the dilaton and the moduli, and at the 
same time provide a mechanism of supersymmetry breaking. The favorite mechanism of such a type is due to 
gaugino condensation in the gauge group hidden sector (see, for instance, p5| , |8j| , and references therein). 
The resulting potential for the dilaton is a combination of exponentials and polynomials in the field. A detailed 
investigation of these condensate models has demonstrated the need for at least two condensates to form if the 
dilaton potential is to develop a minimum at a realistic value ("racetrack" models) [jB6| . Recently an alternative 
proposal, relying on only one gaugino condensate, has been suggested as a method of obtaining a minimum for the 
dilaton |8q| , [jS9f . In this scenario the Kahler potential requires string desired nonperturbative corrections. The 
analysis of this models indicates that it is possible to have a minimum with zero or small positive cosmological 
constant |]90f . The possibilities for inflation and attempts for the resolution of cosmological moduli problem in 
the context of these mechanisms have been considered in J90|-|)3) (and references therein). 

Another possibility of a relaxation mechanism by which various moduli fields are attracted towards their 
present vacuum expectation values due to string-loop modifications of the dilaton couplings offers Damour- 
Polyakov mechanism |29|]. The idea is that nonperturbative effects associated with higher genus corrections, 
may naturally generate different nonmonotonic coupling functions of the dilaton to the other fields (about the 
possibility of stabilizing the moduli fields of type II string theory or M-theory by using the tree-level couplings to 
the non-trivial NS-NS and R-R form fields see |32|). Under the assumption that the different coupling functions 
have extrema at some common point, the expansion of the Universe drives the dilaton vacuum expectation value 
towards the extrema at which the interaction with matter become very weak. Constraints for this mechanism 
from the big bang nucleosynthesis and from the emission of massless dilatons by the binary pulsar are studied 
in |54| . The authors of |53| have shown that the inflationary era with Damour-Polyakov mechanism could solve 
the cosmological moduli problem and the produced quantum fluctuations of the relevant moduli fields during 
this era are naturally compatible with the observational requirements. 

In |66| for the general case of the dilaton coupling functions we have derived the solutions to the string 
cosmology with arbitrary curvature and with moduli and Kalb-Ramond fields as a source. These solutions 
were given in terms of integrations depending on the coupling functions. For the general case of the barotropic 
perfect fluid as a non-gravitational source the set of cosmological equations can not be solved explicitly. Here 
we shall present the phase-space analysis of the string cosmologies with loop-corrections to the dilaton coupling 
functions and barotropic perfect fluid. The issue of dilaton stabilization is considered within the framework of 
Damour-Polyakov mechanism. 

We have organized the paper as follows. In the next section the string effective action with dilaton coupling 
functions of general form and its different conformal frames are considered. Section 3 concentrates on homo- 
geneous and isotropic cosmological models. For the barotropic perfect fluid the main features of the model are 
determined by the barotropic index and by the function a(<p). At tree-level when the dilaton coupling functions 
are exponentials for the NS-NS sector fields and constants for the R-R sector this function is a constant. In 
the general case of the models with curved space the set of cosmological equations can be presented in the 
form of third-order autonomous dynamical system. In section 4 the cases are considered when this system 
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is exactly integrable. They include pure gravi-dilaton, stiff fluid, radiation-dominated models, as well as two 
tree-level special solutions for the general barotropic index. The section 5 presents the qualitative analysis for 
the spatially-flat models, when the corresponding dynamical system is reduced to the second-order one. For 
this system the critical points and their stability are investigated. Various qualitatively different phase portraits 
are presented. Another case when the corresponding dynamical system is reduced to the second-order system 
are tree- level models (section 6). The spatially curved models in the general case of dilaton coupling functions 
are considered in section 7. For these models the phase space is three dimensional. All possible critical points 
are found and their stability is investigated. The phase portraits are plotted for both closed and open mod- 
els. The qualitative evolution of general solution at infinite values of the phase space variables is investigated. 
The dilaton stabilization within the framework of Damour-Polyakov mechanism is considered. Our results are 
summarized in section 8. 



2 String effective action 



The beta functions for strings propagating in a background of massless fields are the equations of motion of a 
certain master spacetime action which can be computed as an expansion in the string tension a' . The general 
form of the corresponding D - dimensional effective action, including possible loop corrections, can be written 
as H-S §> 1, H 



S 



G 



-F R (tp) R - AF V (tp) d M tpd M tp + L m ltp, G MN , $ 



(2.1) 



where tilted letters specify the quantities in the string conformal frame, tp - is the dilaton field, R denotes 
curvature scalar of the D - dimensional metric Gmn, and ip stands for all other degrees of freedom. Including 
the most important special cases, the Lagrangian density L m has the form 



V F A (<p) 
4-2(5 r + l)! 



F 



Y,F Xi (tp) G MN d MXi d NX i 



(2.2) 



Ffr (tp) Ipi Dlpi - TTliFmt (tp) Ipi tpi 



Here F r — dA r ,r — 1,2, ... are the strengths of the 5 r - forms A r of both R-R and NS-NS types, Xi are scalar 
fields (including various moduli associated with compactification of extra dimensions), tpi stand for fermion 
fields. The potential term U includes, in particular, a possible nonperturbative dilaton potential. Note that 
in the supergravity theories deriving from superstring theories one naturally encounters form fields of ranks 
1-4. For example, the 101? type IIA string action contains the universal Kalb-Ramond 2- form field Bmn with 
strength Hmnp in the NS-NS sector, and 3-form field coming from 111? supergravity and second rank field 
strength for the Kaluza-Klcin vector Am in the R-R sector. 

At the present stage of development of string theory we do not know the global behaviour of the dilaton 
coupling functions Fr (<p) beyond the fact that in the weak coupling limit (tp — ► — oo) they should admit an 
expansion 



F K (tp) = e 



-2CV 



1 



£4^ 

i=i j 



(2.3) 



The first term in this expansion is the string tree-level contribution. At this level the dilaton coupling is a 
uniform e~ 2ip (£ = 1) in the NS — NS sector, and it does not couple (in string frame) to the R — R sector field 
strengths ( £ = ). The dimensionless coefficient Z$ represents the I -loop contribution, and the parameter 
of the loop expansion is e 2v . Note that in the low-energy regime, with broken supersymmetry there are no a 



priori obstacles to having couplings (2.3) with Z K ^ 0. 

The effective action (2.1) is written in the form directly obtained from string a - model calculations (so-called 
string frame). One may rewrite it in a different frame by making dilaton depended rescaling field redefinition 



G 



MN 



= n- 2 (ip)G 



MN 



(2.4) 
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By choosing the scale factor as 



n = q e , n E {(p) = Fj/ (1 - n) 



(2.5) 



one obtains the Einstein (E-) frame action, where the graviton kinetic term is diagonalized. The corresponding 
action takes the canonical form 



S = 



J d D xy/\G\ [-R - AF V (<p) d M <pd M <p + L m (ip, G MN ,iP)] , 



(2.6) 



where the indices are now raised and lowered with Gmn- The function in front of dilaton kinetic term and 
nongravitational Lagrangian density are related to the corresponding string frame functions by 



4^ (ip) 



F in 



1 \F, 



VL D L,, 



(ip,Vl 2 GMN,lf) 



(2.7) 



For Ftp > the dilaton field would be ghost-like, since its kinetic term would be negative. For this reason we 
shall consider the case F v < . Under this assumption by introducing a new scalar field <fi as 



d<j> — 1\J '—Ftpdtp 

the dilaton kinetic term can be written in canonical form. At tree-level F v 
field is proportional to the dilaton: 



\AT- 1 

In the general conformal frame related to the string one by transformation (|] 
also can be written in canonical form rescaling the fermion fields (see also pjJ ) 



The corresponding Lagrangian takes the form 



1/2 



L f = 



E 



TrHil(y>)-~ — —ViVi 
F^ W) 



(2.8) 

T/(n — 1) and the new scalar 

(2.9) 

|), the kinetic term of fermions 
(2.10) 

(2.11) 



In the E-frame we have to substitute the conformal factor from ((2^ 

Another important conformal frame of scalar-tensor theories is the so-called Jordan frame, in which the 
nongravitational part of the action does not contain the scalar field. In this frame the laws of evolution of the 
nongravitational fields take the same form as in multidimensional GR. For example, in the case of fermion field 
ipi , as it follows from ( [2.11 ), the Jordan frame corresponds to the choice 



n = n[P = f^ (if ) /F mi [<p) 



(2.12) 



This frame is very natural, because fermion matter couples directly only to this metric, particles have constant 
masses and move on geodesies of its metric. As we see from (2.12) in the case of Damour-Polyakov universality 
ansatz for the dilaton coupling functions, for fermions the string and Jordan frames coincide. 

For the action (2.1) the Jordan frame, in general, can not be realized. Let us consider the important special 
case, when in (2.1) the dilaton dependence is factorized: 



L, 



Up, G M N,i^j = F L (ip) L (g M n, ip^j ■ 



For these Lagrangians the Jordan frame exists, if the function L has a certain conformal weight (5: 

L(Q?G MN ,ip) =Q^L(G M nA) 



(2.13) 



(2.14) 
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In general conformal frame the corresponding function is determined from the last relation of and has the 
form 

L m = n D+2 ^F L (cp) L (Gmn, i>) (2.15) 
Now the choice of the conformal factor in accordance with 

-l 



fl = fl 



D+2P 
J 



(2.16) 



leads to the Jordan frame with exception of the case f3 = —D/2 , when L m is conformal invariant. For the 
r-forms considered above the parameter /3 = —r — 1 has the following values: [3 = —3 for Kalb - Ramond field, 
(3 — — 2 for gauge field. The scalar part of the Lagrangian (2.2) has a certain conformal weight {j3 = —1) if the 
potential terms are absent and the condition (2.13) is fulfilled if the functions at kinetic terms are universal. 
Assuming the universality ansatz for the dilaton coupling functions we see from (2.5) and (2.16) that the Jordan 
and E-frames coincide for massless scalar fields. In particular this is the case at tree-level. Note that for the 
R-R fields the string and Jordan frames coincide at tree-level. 



3 Cosmological model 

In this paper we shall consider homogeneous and isotropic cosmologies. The corresponding E-frame metric has 
the form 

ds 2 = N 2 (t)dt 2 - R 2 (t)ds 2 n (3.1) 

where ds n is the line element of a n = D — 1 - dimensional space of constant curvature, N(t) and R(t) are the 
lapse function and scale factor. The choice of N(t) corresponds to the various time coordinates. For example, we 
have N = 1 and N = R(t) in the cases of synchronous and conformal time coordinates, correspondingly. From 
the homogeneity of the model it follows that the dilaton field should also depend on time only, ip = <p(t). From 
the field equations we obtain that the energy- momentum tensor corresponding to the metric (3.1) is diagonal 
and can be presented in the perfect fluid form 

T^ = diag(s,...,-p,...) (3.2) 

where s is the energy density and p is the effective pressure. If the corresponding nongravitational Lagrangian 
does not depend on derivatives of the metric tensor, the values of these quantities in various conformal frames 
are related by 

T» = n D T^ (3.3) 
with Q = £Ie,J for the E- and Jordan frames respectively. The functions N and R are frame dependent too: 

N(t) = rr 1 ^)^), R(t) = fr 1 ^)^*) (3.4) 

The E-frame evolution equations for the scale factor and dilaton field can be written as 

H+H(nH-N/N^ = N 2 be-N 2 k(n-1)/R 2 (3.5) 
'(j> + (t> (nH- N /N) = N 2 ae/2 

2 -n(n - l)H 2 = -N 2 e + N 2 kn{n-l)/R 2 

where the overdots denote time derivatives, k is the sign of the curvature of the spatial sections (k — — 1, 0, 1 
for spatially- flat, closed and open models respectively) and the following notation are introduced 



1 SL rn ^\G\ 

a = p/e, a= — = — (3.7) 

ey \G\ 0(f) 
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As is well known, the last equation of ( |3.5| ) is the constraint equation, which is a consequ enc e of first two ones 
and gives only a restriction on the initial values. As we see from the second equation of (3J5) for the solutions 
with constant dilaton = O , the values 4>o have to be roots of the equation a(<fi) — 0. These solutions coincide 
with those of GR. The existence of the function a((f>) with appropriate zeroes provides a mechanism stabilizing 
the dilaton. So far several schemes have been suggested to generate such a function. The most popular among 
those is the generation of dilaton potential due to non-perturbative eff ects , such as gaugino condensation in the 
hidden sector. The corresponding cosmological equations follow from (3.5) by choosing L m = — V(ip). Now the 
E- frame functions e and a in RHS of equations (3.5) are defined as 



V(<p)= F R (<p) V{<p) 



ae 



-V(<p) 



(3.8) 



The corresponding qualitative analysis for spatially flat models have been carried out in |57| ] without specifying 
the form of dilaton potential. 

An alternative mechanism stabilizing the dilaton at phenomcnologically acceptable values have been sug- 



gested in 29 1 . In this mechanism (referred as Damour-Polyakov mechanism) the corresponding function a(<f>) 
capable to stabilize the dilaton, is generated due to the higher-loop corrections (see (2.3) to the functions of 
dilaton coupling with the matter sector. Below we shall consider just this case. 

As a consequence of the dilaton dependence of Lagrangian L m the energy-momentum tensor is acted upon 
by a dilaton gradient force: 

-eadMfj) (3-9) 



rpN 

1 M:N 



In the cosmological context this equation takes the form 



£ je + n(l + a)H + a 4>= 



(3.10) 



For the Jordan frame a = and energy-momentum tensor is covariantly conserved with respect to its metric. 

Note that, in general, the quantities a and a are functions of time. The function a is frame independent, 
and the values of the function a in the different conformal frames are related by 



«s) = a — (1 — na)d In il/dcf>, 



(3.11) 



where an - corresponds to the conformal frame, defined by the metric (2.4). In the Jordan frame the nongrav 



itational Lagrangian does not depend on dilaton field and therefore in this frame a = aj = 0. It follows from 
here that in the case of sources for which the Jordan frame can be realized, in the string and E-frames we have 
respectively 

a = (l-na)dlntlj/d(f>, a = (1 - na) din (Oj/Qe) /# (3.12) 

Note that the c oefficient in these formulas, 1 — na = T/e is equal to zero for the radiation. At tree-level the 
functions ( 3.12 ) are constant and in the E- frame 



1 — na 



1 



(n-l)C 
i + 2f3+l 



- 1 



(3.13) 



1, C = 1) horn (3.13) one has a — and at this approximation Jordan 



For NS-NS massless scalar fields (/? 
and E-frames coincide. 

For the field strength F r in ( |2.2| ) there are two types of ansatze which are compatible with the symmetries 
of the metric (3.1). For the first one, known as elementary, one has r = n and 



FoM-L-Mn = const ■ NR n F A 1 (ip)e Ml ---M„ 



(3.14) 



This form of antisymmetric tensor directly follows from the corresponding field equations. In this case the 
Bianchi identity is vacuous. In n = 3 such an ansatz can be realized by 4-form field strength coming from 11D 
supergravity. Note that for this field the additional Chcrn-Simons term is contained in the Lagrangian (2.2). 
But for the cosmological backgrounds we consider this term is always zero. 



G 



The second ansatz is known as solitonic and is realized by forms with r = n — 1 and 



Fmi-m u = const ■ £mi---m„ 



(3.15) 



For the case n = 3 as an example we can consider the Kalb-Ramond field. In (3.14) and (3.15) the symbol 
£Mi---M„ takes the values or 1 and is completely antisymmetric on all indices. By calculating the corresponding 
energy-momentum tensors it can be easily seen that in the case of solitonic ansatz the form held is equivalent 
to the perfect fluid with equation of state p = e (a = 1), and the form held with elementary ansatz corresponds 
to the perfect fluid with p = — e (a = —1). By taking into account that (3 = — n and (3 = — (n + 1) for solitonic 



and elementary ansatze, from ( 3.13 ) one finds 



[(n - 1)C + n - 



(3.16) 



where e = 1,-1 for the solitonic and elementary cases correspondingly. Note that the tree-level cosmological 
backgrounds with various types of NS-NS and R-R form fields were considered in 0, J27), [fU, fH> 

M, l6l, [fzl, 



In the paper 1 66 for the general case of dilaton coupling functions we have derived the cosmological solutions 
with arbitrary curvature and dilaton, moduli and for the Kalb-Ramond fields. These solutions were given in 
closed form in terms of integrations depending on coupling functions. In present paper by using the qualitative 
methods we shall investigate the main properties of cosmological models without the specifying the dilaton 
coupling functions and for the general case of the nongravitational source with barotropic equation of state. In 
particular, the possibility of dilaton stabilization will be considered by using the Damour-Polyakov mechanism. 



4 Dynamical system and special solutions 



In the following it will be convenient to work in the conformal gauge, when N(t) = R(t). The corresponding 
time coordinate we will denote by index c. By introducing a new independent variable r and the function h(r) 
in accordance with 

/T7^ r„ , dlnR H TT d\nR 

aV = V H 2 + kdt c , h= — — = — — — — H = — - — (4.1) 



dr 



dt c 



(from the last equation of (3J5) it follows, that for the models with e > the expression H 2 + k is nonnegative) 
and substituting N 2 e from the last equation into the first two ones, the equations of motion ( |3.5|) can be 
presented in the form of third-order autonomous dynamical system 



dr 
dh 



x, — — = n(n — 1) 
dr 

(1 - h 2 ) [(n- l)(nb- 
This system is invariant under the transformations 

t — > — r, x — > —x, 
relating expanding and contracting models, and 



-x 2 } [a(<f>)/2-bhx], 
1) - bx 2 ] 



(4.2) 



(4.3) 



(4.4) 



relating the models, for two, in general different, functions a(<p) and —a(—<j>). For the source, satisfying the 
conditions (2.14) and (2.15) with constant a the function a is determined by the relation (3.12) and depends 
only on scalar field <p. In the following we shall consider this case. The energy density can be expressed via the 
solutions of the system (4.2) as 



R 2 e 



(H 2 +k) [n(n- 1) ~x 2 



(4.5) 
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It follows from here that the phase trajectories describing the models with nonnegative energy density lie in 
region 

n(n - 1) - x 2 > 0. (4.6) 

If t = J Rdt c denotes the E-frame comoving time, then for the corresponding acceleration from the first of 
equations ( |3.5| ) one obtains 

'' 2 U TJ2 " k [{n-l){nb-l)-bx 2 ] (4.7) 



dt 2 



R 



As it follows from here when 



a < 2/n — 1 

one has accelerated expansion (in the E-frame), d 2 R/dt 2 > , in the region 

bx 2 <l- n(l +a)/2 



(4.8) 
(4.9) 



The segments of phase trajectories lying in this region describe an inflationary type evolution. 

In the condition ( |4.6| ) the sign of equality corresponds to the pure gravi-dilaton models (e = 0). In this 
case the set of cosmological equations is exactly integrable in terms of conformal time and the corresponding 
solution have the form Ol, [fill, 



R — R n 



Vk(n- l)t c ] /Vk 



l/(n-l) 



±ln 



tg 



Vk(r 



l)t r 



/Vk 



n- 1 



—F v (ip)dip 



(4.10) 



(4.11) 



Given the dilaton coupling functions the second of these formula determine the dilaton as a function of the 
conformal time. E-frame synchronous time coordinate is defined by t — J Rdt c . This integral can be expressed 
through the elementary functions for the spatially-flat models only and the corresponding solution 



R~\t\ 



l/n 



const ± yl — l/n In \t\ 



lies in the basis of pre-big bang inflationary model proposed in 15 
parameter r is related to t by 



(4.12) 

Note that for this solution the 

(4.13) 



The equation of pure gravi-dilaton phase trajectories of the dynamical system (4.2) has the form 

x = ±\J n(n — 1), h — =p tanh fe yjl — l/n ((f) — <pi) 



(4.14) 



with <j)\ being an integration constant. These models represent an invariant two-dimensional subspace of the 
three-dimensional phase space of (4.2). In figure 1 the corresponding phase portrait for the models with increas- 
ing dilaton (upper sign) is plotted on the phase plane (</>, h) compactificd onto rectangle (y, w) in accordance 
with 

0<y<l (4.15) 



i-v' 



w 



H 



1- w 



u= ; -l<to<l (4.16) 

\H\ + VH 2 + k ~ ~ 

Note that dw/dh = (1 + |^|) -2 and the function w = w(h) is monotonic. The horizontal segments w = ±1/2 
(h = ±1) correspond to the E-frame expanding and contracting spatially-flat models. They split the phase 
space into three invariant subspaces corresponding to the closed (\w\ < 1/2), expanding (1/2 < w < 1) and 
contracting (—1 < w < —1/2) open models. The segment h = — 1 represents the pre-big bang inflationary 
solution. At tree-level this solution in the string frame describes an expansion of superinflationary type (pole 
inflation). As we see from figure 1 it is the late-time attractor for k = ±1 contracting (in the E-frame) gravi- 
dilaton solutions. For these solutions inflation in the string frame starts when trajectories approach the segment 



h = — 1 sufficiently close and they have finite time to inflate. To inflate sufficiently to solve horizon and flatness 
problems of standard cosmology the solutions with k = ±1 have to approach the segment h = — 1 at extremely 
early times. This illustrates the fine tuning problem of the pre-big bang models considered in J70[,||73[. 

For the sources with the equation of state p = e (for example, massless scalar field) one has b = and 
the first of the equations (3.E) coincides with the corresponding vacuum equation. As a consequence for these 
sources the E- frame scale factor is described by the same expression ( 4.10| ) (the corresponding solutions for the 
dilaton see in [ |66| ). Another important case when the dynamical system ( |4.2| ) admits exact solution are the 
radiation-dominated models. 



4.1 Radiation-dominated solution 



For the radiation a = 1/n and therefore according to (3.11) a — 0. Integrating two last equations one has 



x 2 = 



n(n — 1) 



1 



h 2 = 1 - 



kuou 2 
1 + u 



Ro 



n-l 



(4.17) 



where uo, Ro are integration constants. As we see for the models with k = 1 the possible values of the function 
u are limited: u < u mi where the maximum value can be expressed via constant uo by means of relation 
u m — (l + \A + 4uq) /2uq- The expressions for the Hubble function and dilaton have the form 



\H\ = 



Vl + u — kuou 2 



4> = 2 J \J ' —F ip (<p)d<p — const ± ^ In ( 



1 + - + 2i?V^o 
u 



(4.18) 



(4.19) 



For the expanding models the upper (lower) sign in (4.19) corresponds to 4>< (> 0). In the case of spatially 
flat models (fc = 0) and n — 3 the solution ( |4.19 ) for the scalar field coincides with that given in [101] for the 
scalar-tensor theories. Note that as u — » oo (which is possible for open and flat models only) <j> move only a 
finite amount and H — > 1 for open models and H — > for spatially flat ones. The scale factor dependence on 
the conformal time is determined from 



2uou 




y/k(n - 1) (t c - 1 ) 
(n - 1) (t c - to)} , 



/Vk, 
k = 



u > -fc/4 
■l,u < 1/4 



.20) 



(n - l) 2 (t c - to) 2 /(4uo)-l 







(when uo = 1/4, k = — 1 one has the solution u = — 2 + 2exp [(n — 1) (t c — to)] ). The function 4>(t c ) can be 
found either by inserting (4.2C) into (4.1£) or more simply by direct integration of <p = cons t ■ J dt c /u with the 



help of ( 4.20| ). The last relation follows from the first integral of the second equation of (3.5) in the case of 
radiation, <f> = const. The E-frame synchronous time coordinate t can be expressed via u as 



t = 



u Ro 



uV(«-l)d u 



kuou 2 



(4.21) 



Note that for closed models the segments x = 0, \h\ < 1 and for open models the rays x — 0, \h\ < 1 are solutions 
of the system (4.2) and represent the corresponding GR solutions. If at the beginning of the radiation-dominated 
era u — u r , then the total variation of the dilaton field in this stage is equal to 



= ± 



In- 



2\/r 



u uf: 



u r (l + 2 v / u^) 



(4.22) 



As it follows from here if, at the beginning of the ra diati on era the value of x is not too close to \Jn(n — 1), the 
total shift of the dilaton is of order of one (see also 101]). 
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The scale factor and comoving time coordinate for the corresponding radiation-dominated solution in the 
string frame can be found by conformal transformation and have the form 



R — Rq 



l/(n-l) 



t = ±- 



R(u)du 



Vl + U — UqU 2 



(4.23) 



For the given coupling functions F^{ip) and F^(if) in the action (2.1) we can obtain the function ip(u) by 
inverting the equation (4. IE) and then derive scale factor R(u) and comoving time t(u) from (4.23). Therefore 
the formulae ( 4.19| ) and (4.23) determine the string frame radiation-dominated solution in parametric form as 
a function of u, provided we specify the dilaton coupling functions. By choosing F v — Fr — e~ 2tp one obtains 
the corresponding tree-level solutions. 

4.2 Special solutions for a = const 



In the case a = const the dynamical system (4.2) admits two classes of special solutions. They are important 
on account of that they are early or late-time attractors for general solution, as we shall see below. For the first 
of these solutions 

x = a/2bh, h = ±l (4.24) 

It corresponds to the spatially flat models and has positive energy density if |a/26| < yjn{n — 1). The depen- 
dence of the scale factor and scalar field on the E-frame comoving time is given by expressions 



R = const ■ \t — to 



ib/cti 



2a 



n-1 



(1-a 2 ) 



(4.25) 



<p = const H In It — to I 

a>i 

Here t > to (t < to) corresponds to the expanding (contracting) models. This solution describes an extended 
inflation when 

a 2 /2b< 2-n(l + a), (4.26) 



the necessary condition of which is (4.8). The corresponding energy density has the form 

,2 



4 a 2 



a\ (t - t ) : 



r, ao = 2byJ n(n - 1) 



(4.27) 



The particular case of models with a = const are trec-lcvcl models, for which the value of a is determined 
by (3.13). In this case the corresponding string frame solution can be obtained by conformal transformation 



R 



i^^^R, t = / e^'^^dt 



and has the form 



R = const ■ \t — to 
= const + 



(4.28) 
(4.29) 



ay/n — 1 



a + \/n — lai/2 



In i — £q 



when \Jn — la ^ — 1 ± \J\ — n(l — a 2 ), and the form 



R = Ro exp 



(l + 2bVn~ I/a) Ht 



V = -(n-l)H(t-t Q ) 



(4.30) 



otherwise (note that in the spatially-flat case the tree-level models are exactly integrable in parametric form 
for general anisotropic barotropic fluid (see |f47jl)). Here t is the string frame comoving time, related to the 
parameter r and E-frame comoving time by 



r/(4bh) =ar 1 ln|t-t | = [2aVn - 1 + m] In |t — t | 



(4.31) 
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For the considered special solution the relative contribution of dilaton to the energy density is equal to a 2 /oq. 
When tree-level dilaton coupling to the nongravitational matter is trivial in the string frame (the case C = in 
( |2.3| )) from (4.29) one obtains the solution previously considered in [p4| . The analogous solutions in the case 
C = 1 are investigated in pTJ . Note that for the case 



1 -a 



[2-n(l + o)] 



(4.32) 



the scale factor is linear function on comoving time in both E- and string frames. As it follows from (4.29) when 



1 - a 



(4.33) 



one has solution with constant string frame scale factor and logarithmic or linear dilaton in the cases ( 4.29p 
and (4.30) respectively. For example, this condition is satisfied in the pure gravi-dilaton model of noncritical 
strings, when in the action (2.1) L m = (D — D cr ) e~ 2ip /3 , D cr - critical dimension ( =10 for superstrings). In 
this case a = — 1,(3 — and in accordance with ( 3.1 3| ) a = —2/y/n — 1 , which coincides with (4.33). This 
solution corresponds to the one of (4.30) with linear dilaton and was considered previously in Jl2| . 
For the second class of special solutions of ( |4.2[ ) with a — const one has 

(4.34) 



x = ±\j(n - l)(n - 1/6), h = a/2bx 



The necessary condition for existence of this solution is (4.8). In the E- frame the corresponding time dependence 
is determined by the expressions 



R = H(t-t ), H = ± 



n(l + o) . 
2-(l-a) 



(n — l)a 
const + a _1 [2 — n(l + a)] lni? 



-1/2 



(4.35) 



with to being integration constant. Here upper/lower sign corresponds to the expanding/contracing models. At 
tree-level the same solution in the string frame can be obtained by using the conformal transformation ( 4.28| ) 
and takes the form 



R = H(t- t ), cj> = const + [1/y/n- T + a/(2 - n(l + a))]^ 1 lni? 



In the case 



a — [n (1 + a) - 2]/Vn- 1 
one has a solution with constant string frame scale factor instead of ( 4.36| ): 



R = const, 



= ± 



t-t 
i — — — 
R 



(4.36) 



(4.37) 



(4.38) 



In the limit \h\ — > 1 in ( 4.34 ) this solution tends to the previ ous special solutions ( 4.24 ) corresponding to the 

spatially flat models^ Note that for the values a from ( 4.37 ) one has h\ = 1 — (n — 1)(1 + a)/(l — a), and 

hence the solution ( 4.38 ) corresponds to the c losed models when a > — 1 (recall that f or the solution ( 4.34 ) the 
condition (fh^) must be satisfied). If in ( [4.37] ) a = —1 (a = —2/y/n — 1) the solution ( 4.34 ) corresponds t o the 
spatially-flat models (h — ±1) and has the form ( 4.30]) with the zero exponent. Note that the solution (4.34) 
corresponds to closed models when the condition (4.26) is satisfied. As it have been mentioned above in this 
case the solution (4.24) describes an extended inflation. 
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5 Qualitative analysis of spatially-flat model 



For general case of the parameter a the system ( |4.2| ) cannot be solved explicitly. To obtain the generic features 
of the cosmological evolution we shall use the dynamical systems methods. In this section we shall focus on 
generic cosmological properties for the spatially-flat models. For the cosmological model with flat space (k = 0) 
from (4.1) one has h = 1,-1 for the expanding and contracting models respectively. Now ( |4.2[ ) is reduced to 
the second-order dynamical system 



d4> 

rf7 



— = \n(n 
dr L 



1) - x 2 ] (a/2 - bhx) 



(5.1) 



const -hlriR, h = ±1 



Below we shall consider expanding models. The corresponding results for the contracting ones can be obtained by 
transformation t — ► — r, x — ► —x. The trajectories of the system (5.1), describing the solutions with nonnegative 
energy density lie in the horizontal strip (|4.6|) of the phase plane ((f), x). The boundary of this strip 



± v / n(n-l) 



(5.2) 



is a solution of the system and corr esponds to the pur e gr avi-dilaton solutions (4.10), (4.11). At tree-level the 
function a is a co nstan t (see ( 3.13|) ) and th e sys tem (5.1) has another special solution ( 4.24 ) with comoving 
time dependence ( 4.25 ) in the E- frame and ( 4.29|) in the string frame. 

For the general tree-level spatially flat solutions the equation for the phase trajectories have the form 



a 



= const + In 



a 
2b 



\J n(n — 1) + . 



\J n(n — 1) — x 



(5.3) 



±oo the general solution tends to one of 



where «o is defined b y ([4.271). A s we see in the limiting regions i 
the special solutions ( |5.2|) , ( [4.24| ). 

For the case of radiation as a nongravitational source we have a — l/n and in accordance with ( |3.12 ) a = 0. 
In this case the system ([Tl]) has a solution with constant dilaton, which coincides with the corresponding GR 
solution. The equation of the phase trajectories describing the radiation-dominated solutions with variable 

dilaton has the form 

x = -^Jn(n - 1) tanh [a<y((j) - 0i)/2] , (5.4) 

with 0i being an integration constant. The corresponding r dependence is given by the expression (the analogous 
solutions in scalar-tensor theories of gravity for n = 3 see 101 ) 



= 01 ±2 



1 



In 



e _(n-l)(T-T )/2 + v / 1 + e -(n-l)(r-T ) 



(5.5) 



with a new integration constant r . 

We now turn to the systematic analysis of the qualitative structure of the phase portraits of dynamical 



system (5.1) for the general case of the function a(<j)). The key aspect for this is to determine the critical points 
and their stability. For the system (5.1) those are the following points of the phase plane (<p,x): 



(00,0), a(0 o ) = O 



(5.6) 



Note that according to the relation (3.12) the zeroes </>o are extremums of the function flj/flE- Here Qe have the 
form (2.5) and for the sources with the Lagrangians of given conformal weight the function £lj is determined 
by (2.16). In particular by assuming the Damour-Polyakov universality ansatz for the higher genus dilaton 



couplings, one finds that the zeroes (5.6) are extremums of the dilaton coupling function. 

Near the critical points the corresponding linear system coincides with the equation of damped oscillator. 
The character of the critical points is defined by the eigenvalues 



Ai, 2 =-k ±yfkJ+ n(n - l)a /2, a' = (^j 



k 



-(1-a) 



(5.7) 



For the real Ai the angula r coe fficients of the separatrixes in this point are equal to ki — \i. The standard 
analysis (see, for example, |l03| ) leads to the following results. 
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1. When a' > the eigenvalues Ai are real and have different signs. The critical point is a saddle (see figure 
2f). 



In the case a' < the critical point (5.6) is stable (asymptotically stable for a > 0) for the expanding 



models and unstable for contracting ones. It is necessary to distinguish the following cases. 

2. When ao = , a' Q < the critical point is a center (see figures 3e and 3f). Near this point the solution 
for the scalar field has the form 

4> — 4>q = 4>\ cos (loqt + 9) , =—n(n— l)a' /2 (5.8) 

This case corresponds to sources with the equation of state p = e, as it is the case for massless scalars 
and form-fields with solitonic ansatz. For D = 4 the antisymmetric Kalb-Ramond field is reduced to the 
pseudoscalar axion field as well. The corresponding solutions are considered in p6[ . 

3. a' Q < — <Xq/8. The critical point is a stable sink (see figures 3a and 3d) and the solution in its vicinity 
represents damped oscillations for the expanding models: 



ie - fc0T cos + (5.9) 



4. a' = — Q!q/8. The eigenvalues are real and equal to —^Jn(n — l)ao/4. The system is critically damped 
and the corresponding point is an improper node (see figure 3b). Near of this point the solution has the 
form 

(j) - 0o = (oi + a 2 r)e- kT (5.10) 

with a,i being real constants. 

5. — «o/8 < a' Q < 0. The eigenvalues are real and negative. The critical point is a stable node (figure 3c). 
As | Ai| < | As | near of this point all trajectories except special ones, touch the eigenvector corresponding 
to the eigenvalue Ai and the solution looks like 

4>-(/)o = b ie XlT + b 2 e X2T (5.11) 

6. And at last let us consider the degenerate case a' = 0, when Ai = 0. Let near of the critical point the 
function a(<p) has an expansion 

= A«o m) (0 - M m + ■ ■ ■ , m>2. (5.12) 
m! 



The standard analysis (see, e.g. 102]) leads to the following results 



(a) ttl = 2. The critical point is an equilibrium state with a node sector (stable for h = 1, and unstable 
for h = — 1 ) and with two saddle sectors. The angular coefficients of the separatrixes are equal to 
k% = 0, &2 = —\Jn(n — l)ao/2. In the case a ' > 0(< 0) the trajectories of the node sector tend to 
the critical point from the left (right) of the separatrix with coefficient k 2 - 

(b) Tji = 3. The point has a character of a saddle (node) when 

1 + 4a 3) /3ao > 0(< 0) (5.13) 

(c) m > 4. The critical point is a saddle. 



Uniting the expressions (5.8), (5.9), (5.11) in the general formula 

<t>- 4> = 4> 1 e aiT cos{a 2 r + a 3 ) (5.14) 
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with appropriately defined constants aj, from the first of the equations (3.5) it is easy to find the behaviour of 
the Hubble function near the critical point 



H = 



n(l + a)t 



{1 + const ■ e 2aiT [ci + c 2 sin(2a 2 r + c 3 )]} , • 



±2 In jt| 
n(l + a) 



(5.15) 



where upper/lower sign corresponds to the E- frame expanding/contracting models, and constants can be ex- 
pressed via ai and parameters a, n. The first summand in the curly brackets corresponds to the GR solution 
and the second one represents the small correction due to variable dilaton. 

From the above analysis it follows that in the case of expanding models the existence for the function a(4>) 
of zero 4> — <fio with a'(0o) < leads to an efficient mechanism for the dilaton stabilization and relaxation of 
string effective gravity to the GR. Indeed, for the relation of the dilaton variations at some time moments U 
and t„ we have 



4>(U) - 4>o 



exp{- |fleAi| [t(U) - t(U)]} 



R{U) 



R{t e ) 



|i?eAi| 



(5.16) 



In the cases 3) and 4), when critical point is a sink or improper node one has |i?eAi| = n(l — a)/A. In the 
inflationary stage a = — 1 and the relation of the scale factors at the beginning and at the end of inflation is 
~ e -65 . If initially the dilaton shift ~ 1, one receives that at the end of inflation the dilaton variations are 
strongly suppressed (n = 3) (see |5^| ) 

0(t e ) - O < 1CP 42 (5.17) 



Assuming that for the various sources the zero of the function a{4>) is the same (in particular, this is the case 
when the dilaton couplings are universal) the subsequent expansion stages (radiation- and matter- dominated) 
leads to the further suppression of these variations: \5(j>\ < 1CP 49 at present epoch |53|| . Such variations are 
much more far from observational restrictions following from the tests of the equivalence principle. Thus this 
scheme leads to the natural suppression of the dilaton variations in the expanding universe and lies in the basis 
of dilaton stabilizing mechanism known as Damour-Polyakov mechanism. Note that this scheme will work also 
in the case 5), if the value of \a' \ is not too small. 



In order to fully describe the qualitative evolution of the system (5.1), we must investigate the behavior of 
the phase trajectories at infinity of the phase plane. For this it is convenient to map the phase plane onto a 
finite region. This can be done by introducing a new variable y in accordance with (4.15). In terms of this 
variable the system (5.1) takes the form 



dr 



2/(1 - y)x, 



dx 
dr 



[n(n - 1) - x 2 ] (a(y)/2 - 6a;) 



(5.18) 



The phase space (y, x) of this dynamical system corresponding to the models with nonnegative energy density 
is defined by the relations <y <\,\x\ < i/n(n — 1). The critical points (y,x) on the boundary of this region 
are the following points 

(j,±y/n(n-l)) , j = 0,1 (5.19) 



{j,a(j)/2b), \a(j)\<a 
The eigenvalues determining their stability are equal to 



Ai = ±(-l) j y/n(n-l), A 2 = 2bn{n - 1)(1 =F a{j)/ 



a 



in the case of (5.19), and 



A! = (-iya(j)/2b, A 2 = -bn(n - 1) (l - a 2 (j)/a 2 ) 



(5.20) 



(5.21) 



(5.22) 



in the case of (5.2C). From these expressions it follows that the critical points are nodes or saddles. The points 
(0, yjn(n — 1)) and (1, — \Jn{n — 1)) are unstable. Under condition (— iya(j) < the point (5.20) is stable 
node and as 4> — > oo the trajectories tend to the solution described by the separatrix x = a(j)/2b , having the 
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form ( 4.25|) (by replacing a — > a(j)) when t — > oo. In the case a(j) = the eigenvalue Ai = and the critical 
point ( 5.20 ) is degenerate. The additional analysis shows that if near y = j,j = 0, 1 one has 



a(y) 



U)(y-j) m + 



a W (j) =0, i < m 



(5.23) 



(j) < and a 



then for the expanding models the critical point ( |5.20| ) is an stable node when (1 — 2j) m+1 c 
saddle point when (1 - 2j) m+1 a ( . m \j) > 0. 

For the sources with a — 1 (for instance, massless scalars or axions) the only critical points on the phase 
space boundary are the points ( 5.19 ) with eigenvalues ( 5.21 ). In this case b — and A2 = =F-i/n(n — l)a(j). 
When a(j) — all points of the segment y — j, \x\ < \Jn(n — 1) are critical points. 

On the basis of the above analysis one can plot the phase portraits of the dynamical system ([O]) for 
an arbitrary fun ctio n a{<j)). The corresponding phase-diagrams for the function — a(— 4>) can be obtained by 
transformation (4.4). Below in figures some qualitatively different cases of the E- frame expanding models 



(upper sign in (5.1)) are pres ented. The portraits for the contracting models can be obtained from here by 
using the transformation ( [4.3| ). In fi gures 2 and 3 the phase plane (<f>,x) is mapped onto rectangle < y < 1, 
I a; I < y n(n — 1) according to ( 4.15 ). The horizontal segments x ~ ±^n(n — 1) correspond to the pure gravi- 
dilaton solutions (4.10), (4.11). For the point (0, ±-\/n(n — 1)) ((l,±y/n(n — 1))) one has r = —00 (+00) 
and t = +00 (—00) for the upper and lower signs respectively. Therefore in accordance with ( 4.13| ) for the 
corresponding values of the E-frame synchronous time coordinate one has t = 0(+oo) and t — +oo(0). 

The phase portrait of figure 2a corresponds to the tree-level models when a = const and it is assumed that 
< a < ao. Horizo ntal s egment x — a/26 represent the special solution ( 4.24 ) considered in the subsection 4.2. 
In accordance with ( 4.31 ) the initial and final points of the evolution of this solution correspond to the values 
t = to and t — +00 of the E-frame synchronous time coordinate and this solution is the late-time attractor 
for trajectories with e > 0. This special solution splits the phase space into two subspaces. The trajectories 
corresponding to the first one evolve from the weak coupling region (ip <C —1) near the pure gravi-dilaton 
solution x = \/n{n — 1) and finish the evolution in the strong coupling region. For the second class of solutions 
(trajectories lying below of the segment x — a /2b) the evolution starts in the strong coupling region at finite time 
moment of the E-frame synchronous time near the gravi-dilaton solution and finish that in the strong coupling 
region too at t — > +00. When a — > ao the special solution ( 4.24 ) tends and in the limiting case coincides 
with the solution x = v/n(n — 1) . At further increase a > ao the qualitative character of the phase portraits 
corresponding to the models with nonnegative energy density, remains the same. Thus the phase portraits for 
a > ao can be obtained from figure 2a by limitary transition a — > ao and are qualitatively equivalent to the 
phase portrait between the solutions x = a /2b and x = —y/n(n — 1). 

The tree-level phase portrait for the case — ao < a < is plotted in figure 2b. All solutions with x > a/26, 
except the pure gravi-dilaton ones, finish their evolution in the weak coupling region near the special solution 
(4.24), which is the late-time attractor. The solutions with x < (>)a/2b evolves from the strong (weak) coupling 
region near the gravi-dilaton solution with decreasing (increasing) dilaton. Note that now there are solutions 
for which during all evolution dilaton remains in the weak coupling region. When a — > — a the solution 
(4.24) approaches to the gravi-dilaton solution x = —^/n(n — 1). At further decrease a < —ao the qualitative 
character of the phase portrait remains the same and is equivalent to that of figure 2b between the solutions 
x = a /2b and x — y n(n — 1). Note that in the considered cases the pure gravi-dilaton solutions are unstable 
when |a| < ao and are stable when \a\ > ao. 

1 / n) one has a = irrespective to the 
is presented in figure 



From (3.12) it follows that for the radiation dominated models (a 
dilaton coupling functions. The corresponding phase portrait of the dynamical system 
2c. All points of the segment < y < 1, x = are critical and correspond to the solution with constant dilaton. 
For the expanding models the solutions with variable dilaton originate from the unstable nodes (0, \/n(ji — 1)) 
and (1, —y/n(n — 1)) at some finite time moment, and , except the pure gravi-dilaton solutions, tend to the 
corresponding solutions with constant dilaton when t — > +00 . The dilaton limitary value depends on the initial 
conditions. The equation of the phase trajectories is determined by the relation ( |5.3p , and the dependence on 
r - by relation (5.5). 



For the general case of dilaton coupling functions the character of the phase diagrams of the dynamical 
system (5J) depends on the number and stability of the zeroes of a(<p) and on the boundary values a(±oo). 
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When a(±oo) > ao (< — ao), ot(4>) > (< 0) the only critical points are (y — 0,1; x — ±yjn(n — 1)) on the 
phase space boundary. For these cases the phase portraits are qualitatively equivalent to the tree-level ones 
with a — const > ao(< — ao) which are considered above. In figure 2d we have plotted the case a(±oo) < ao, 
a(4>) > 0. The general solution with e > approaches the special solution with late-time asymptotic ( 4.25| ), 



where one has to substitute a — a(y = 1). The case a{— oo) > ao, a(+oo) < ao, ct(4>) > is presented in figure 
2e. Now all models begin and finish their evolution in the strong coupling region and the late-time behaviour is 
the same as in the previous case. When — ao < a(— oo) < and < a(+oo) < a the phase portrait is depicted 
in figure 2f for the function a(</>) with a single zero. At this zero a'(</>o) > and the corresponding critical 
point is a saddle. The separatrixes of this saddle point split the phase space into four qualitatively different 
regions. When a(— oo) < — ao or a(+oo) > ao the corresponding portraits can be obtained from figure 2f by 
the corresponding superposing of the points. 

In all considered examples, apart the radiation dominated case and stable separatrixes on figure 2f, the 
models finish their evolution on the phase space boundary (j> — > ±oo ( y = 0, 1 ). This models cannot stabilize 
field <j) at finite values. In figures 3a-3d we have plotted the phase portraits for some simple models with variable 
a , leading to the dilaton stabilization as a result of cosmological evolution. They are illustrations of Damour- 
Polyakov mechanism. The phase portrait in figure 3a corresponds to the function a(</>) with a(— oo) > ao 
, a(+oo) < — ao and it is assumed that at the zero </> = 4>q of the function a(</>) its derivative is negative 
and satisfies to the condition |a'((/>o)| > — ao/8, when the critical point is a sink (case 3)). For the assumed 
boundary values a(±oo) all critical points on the phase space boundary are saddle points. For all spatially-flat 
models, apart the pure gravi-dilaton ones, the expansion of the universe leads to the dilaton stabilization. When 
< a(— oo) < ap, a(+oo) < — ao the critical point appears on the boundary y — 0, having saddle nature and 



the unstable separatrix of which corresponds to the special solution ( 4.24 ) at <j) — > — oo. The phase portrait of 
this case is plotted in figure 3b, where it is assumed the stable critical point is an improper node, |ao| = — ap/8 
(case 4) ). And at last when < a(— oo) < a , — a < a(+oo) < the new critical saddle point appears on 
the boundary y = 1. The phase portrait is presented in figure 3c assuming that a' Q > — (case 5)), when 
the critical point is a node. In all considered cases for all the models, except special ones, corresponding to the 
separatrixes of critical points ( |5.20 ), the initial stage of evolution is dilaton-dominated and the contribution of 



the nongravitational matter can be neglected. The phase portraits for the other possible combinations of the 



relations between a(±oo) and ao can be obtained from the cases considered above by transformation (4.3). In 
figures 2a-c are presented most favorable (from the point of view of dilaton stabilization) cases, when for the all 
models with e > the cosmological evolution transfers the system into the states with small dilaton variations 
(all solutions tend to the models with constant dilaton). As it follows from the expressions of the eigenvalues 
( p>i"l ), fT22] ) at the phase space infinity there are no stable points when 



a(-oo) > 0, a(+oo) < (5.24) 

These conditions are necessary for that the all homogeneous and isotropic models of effective string theory 
during the cosmological evolution have to tend to the solutions with fixed dilaton. In the case of nongravitational 
sources, for which the Jordan frame can be realized, the boundary value a{— oo) is determined by the expression 
(113). For the potential-dominated source a = — 1,(3 — and hence in accordance with ( 3.13] ) a(— oo) 



-2/y/n — 1 and the point (5.20) with j = is stable node. As in the previous case for the part of solutions 
the existence of the zeroes of the function a((f>) with a' < does not lead to the relaxation toward the GR 
solutions. In figure 3d we present the typical for these cases phase portrait, when — ao < a(— oo) < 0, the 
function a(4>) initially increases with increasing <p, becomes zero with the positive derivative in this point (left 
critical point in figure 3d, being saddle point). Further the function reaches the maximum after which decreases 
to the values — ao < a(+oo) < 0. It is assumed that for the second zero one has |a | > ap/8 and it is a sink. The 
cosmological models leading to the final fixed value of the dilaton begins from the unstable nodes (0, d n(n — 1)), 
(1, \/n(n — 1)) in directions lying above the stable separatrixes of the saddle point, corresponding to the left 
zero of the function a{4>). All the other models evolve toward the weak coupling region and approach to the 



solution (4.24) when t — > +oo. The case a(— oo) < — ao or a(+oo) < — ao can be obtained from figure 3d by 



simple superposing of the points (0, a(— oo)/2b) and (0, —y/n(n — 1)) or (1, a(+oo)/2&) and (1, —y/n(n — 1)) 
respectively. 

For R-R fields when 1 — na > one has a(— oo) < and the point ( |5.20| ) with j = is stable node. Therefore 
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there is a class of the solutions with ip — > — oo when t — > +00 and for these solutions dilaton cannot be stabilized 
at finite values. For massless RR axion a(— 00) = \/n — 1, however in this case, as it mentioned above, the 
zeroes of the function a(4>) with a' < are centers, near of which dilaton oscillates without damping (case 2), 
and therefore these sources can not stabilize the dilaton. This is the case for the sources with the equation 
of state p = e (a = 1) and for those 6 = 0. For the function a{(j)) with a single zero <f) = 4> , a'(4>o) < the 
corresponding phase portrait is similar to that of figure 3a with center instead of sink. In figure 3e we have 
plotted the phase portrait for the more complicated function with two zeroes: a'(0o) > / < for the left/right 
one. For the corresponding limitary values we have a(±oo) < 0. There is a subclass of solutions with periodic 
dilaton. When a — 1 and a(y = j) = 0,j = Q, 1 the all points of the segment y — j are critical. As an example 
in figure 3f we have plotted the case a(y — 0) > ao, a(y = 1) = 0. There are two classes of solutions: with 
periodic dilaton, oscillating near the zero <p = <po, and models which begin and finish their evolution in the 
strong coupling region. 

On the basis of the analysis carried out above it is easy to investigate more complicated cases of the function 
a((f>). In the following sections we shall consider the models with curved space. 



6 Spatially curved models with constant a 



Another important case when ( |4.2| ) can be reduced to the second order dynamical system are models with 
constant a. In particular they include tree-level models. Now the field cf> is not contained into the last two 
equations and we have the following second-order autonomous dynamical system 



dx 

dh 
d7 



[n(n - 1) - x 2 ] (a/2 - bhx) (6.1) 
(l-h 2 ) [(n- l)(nb- 1) - bx 2 ] , 



For the models with nonnegative energy density the phase space (x, h) of this system is determined by the 
condition x 2 < n{n — 1), and h 2 < 1 (> 1) for k = 1(— 1). 

6.1 Critical points 

The critical points of ( |6.l[ ) are the following: 
1. Points {x, h) with coordinates 

(±y/n(n-l), h ) , h = -1, +1 (6.2) 

They represent the pure gravi-dilaton spatially-flat solutions (the horizontal segments in figure 1). The 
corresponding eigenvalues and eigenvectors are determined from the relations 

A2 = \J n(n - Vjho (a =p ^o") , \3 = 2h (n— 1) (6.3) 

n 2 = (l,0), n 3 = (0,l) (6.4) 

As it follows from here the points corresponding to the expanding models ( ho = 1) are always unstable 
and have a character of a saddle point when ao < ±a and character of an unstable node when ao > ±a. 



The directions of separatrixes are determined by the eigenvectors (6.4). In the case of node near of the 
critical point all trajectories, except special ones, touch the special solution h = 1 (eigenvector 112 ) when 
±a > ao — 2y/l — ljn and the solutions x = ±y/n(n — 1) (eigenvector ni ) in the other cases. The 



behaviour of the trajectories near the points ( J6 . 2| ) with ho — —1 is obtained by transformation (4.3). 
2. Points 

(±a/26,±l), |a|<a (6.5) 
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with the eigenvalues and eigenvectors 



A 2 = — n(n — l)bh 1 j 



On 



1 of 

A 3 = -2n(n - l)bh ( 1 - — - — 

7X0 C^n 



n 2 



(1, 0), n 3 - (1, 2b/ a + A/na(a 2 /al - 1)). 



(6.6) 



(6.7) 



They correspond to the special solutions (4.24). The point with upper sign is a saddle when a 2 /a^ > 
1 — 1/nb, the direction of the unstable separatrix of which is determined by vector 113, and stable node in 
the other cases. The character of the point with lower sign can be obtained by taking into account the 
transformations ([Q|). 

3. In the case a < 2/n — 1 one has also the following critical points 

(±v/(n-l)(n-l/fe),±/n(a)) , hi(a) 

with the eigenvalues and eigenvectors 



a y/l - 1/nb 



Xi = ±- 



(6.8) 
(6.9) 

n 2>3 ~ (1, -26A 3 , 2 /a(n - 1)) (6.10) 

and corresponding to the special solutions ( [4.34 ). Note that for closed models h 2 < 1 and as it follows 
from (6.9) the eigenvalues have opposite signs and the point is a saddle. For the open models h 2 > 1 and 
the point (^8) with upper/lower sign and for a > is stable/unstable node (sink) when 



-hx {-1 ± [1 - 8(nb - 1)(1 - K 2 )} 1/2 } , i = 2,3 



8(nb- 1)(1 - h- 2 ) < 1(> 1) 



In the E-frame the point (£ 
In the case 



3) corresponds to the solution ( 4.35 ) considered in the subsection 4.2. 
\a\ = aoyl — 1/nb 



(6.11) 



(6.12) 



in (3.8) we have \hi\ — 1 and the critical point ( |6.8| ) corresponds to the spatially-flat models. Now this point 
coincide with the point ( |6.5| ), and A2 = 0, that is it is degenerate critical point. The standard qualitative 
analysis shows (see, for example, [ 103 1 ) that this point represents an equilibrium state with one stable node 
sector and with two saddle sectors. The separatrix between these sectors is the solution x = ±y/n(n — 1). For 
the node sector h 2 > 1 and the corresponding solutions represent open models (for this case the phase portrait 
of the open models see figure 4f ). The trajectories of the saddle sectors correspond to the closed models, and 
the phase portrait is qualitatively equivalent to the case of that figure 4e. In the case (6.12) the corresponding 
solution have the form 



R = const ■ \t\ , <j) = const + — In \t\ 



(6.13) 



R 



const ■ \t\ , 4> — const + 



ln\t\ 



2b/a+l/Vn^T 



Then a = — 1 we have the solution ( 4.3C ) with zero exponent. Here the solution in the string frame is obtained 
at tree-level approximation, when the dilaton field can be found from the rel ation ip = yj n— 10/2. 
In finite part of the phase plane the only critical points of the system (|6.l|) are points (|6.2|) 



For 



the closed models the phase space of the system (6.1) is compact and there are no other critical points. 

For the models with k = — 1 in addition to the analysis of stability of the critical points carried out above, it 
is necessary to investigate the behaviour of phase trajectories at infinity of the phase space, when h — > ±00. Due 
to the invariance with respect to the transformation (O) it is sufficient to consider the case h > 1, corresponding 



18 



to the expanding models. To compactify the phase space let us introduce a n ew v ariable z = l/h, < 2 < 1. 
In terms of independent variable T, dT = hdr = dhiR the dynamical system (3.1) takes the form 



dx 

IT 

dz 

dT 



\n(n — 1) — a; 2 ] (az/2 — bx) 
z(l-z 2 ) [(n- l)(nb- 1) -6a; 2 ] 



(6.14) 



the phase space of which is the rectangle |x| < \/n(n — 1), < z < 1. Besides the points (|6.2[) , (6.5) and (|6.S|), 
the stability properties of which we have considered above, there are the following critical points of the phase 
space (x, z) 

(6.15) 



[x — ±y/n(n — 1), z = 



with the eigenvectors (6.4) and eigenvalues 

A 2 = 2bn(n - 1), 

and 



A, = 1 



(x = 0, z = 0) 



with eigenvalues and eigenvectors 



A 2 



-bn(n — 
= (1,0), 



1) 



A3 = (n — l)(bn 
2{\/n-a) 



1) 



(6.16) 
(6.17) 

(6.18) 



n 2 = (l,U), n 3 ~ 
The solutions corresponding to this point have the form 

R = t, 4> = const 



a(n — 1) 



(6.19) 



and for them e = 0. By making a coordinate transformation t = (i' 2 — r' 2 ) 1 / 2 , rt = r' it is easy to see that they 
correspond to the flat spacetime. 

The critical poi nts (|6.15 ) are saddles with stable separatrixes x — ±yjn(n — 1) and unstable separatrixes 
z = 0. The point ( |6.17| ) is a stable node for a > 2/n — 1 and a saddle otherwise. For the former case and 
expanding models in a neighborhood of the critical point (3.17) the evolution is curvature dominated and the 
solution have the form 



_2 + 2-n(l+o) 
L cr 

2[l-n(l-a)] 



t 



const 



1 



ac 2 (n-l) f2-n(l+a) 

2-n(l + a) 2(l/ra-a) 



(6.20) 



(6.21) 



with Co and zo being integration constants. Note that the asymptotic formula (6.2C) is the same as in GR. 
When a = 1 — 1/n the solution for the scale fac tor is R w t + c 2 lnt/(4n — 6). In the radiation-dominated case 
the second time depe ndent term in the RHS of ( 6.21 ) is absent and a = 1/n. In the case a — 2/n — 1 (b = 1/n) 
the points ( |6.8| ) and ( 6.17 ) coincide and A2 = 0. An additional analysis of the critical point ( 3.17 ) shows that 
it is a stable node. 



6.2 Qualitative evolution 

Having investigated all possible critical points and their stability we can construct the phase-space diagrams for 
any values of parameters a and a. To p lot ph ase portraits for closed and open models on th e same phase-space 
diagram it is convenient to use w from ( 1.16 ) as a phase space variable instead of z in ( 6.14 ). In terms of w the 
trajectories with nonnegative energy density lie in the rectangle 



\x\ < yjn{n-l), \w\ < 1 



(6.22) 
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of the phase plane (x,w). The boundary segments x — ±y/n(n — 1) correspond to the pure gravi-dilaton 
solutions, and w = ±1/2 correspond to the spatially-flat models. The latter split the phase space into three 
invariant subspaces corresponding to the closed, \ w\ < 1/2, expanding, 1/2 < w < 1, and contracting, —1 < 
w < —1/2, open models. The dynamical system ( |S.14 ) has a two dimensional parameter space (a, a), \a\ < 1. 
In this space qualitatively different regions are separated by curves 



a = 2/n — 1, a = ±ag, a = 0, 



1 



2(n-l) 
rt(l — a) ' 



(6.23) 



where the latter corresponds to the condition hf = 1. 

Below considering various cases of the parameters we shall assume a > 0. The phase portraits for negative 
a can be obtained by taking into account the invariance of the system (5.14) with respect to x — > —x, a — > —a. 
Consider qualitatively different cases in turn. 



1. 



For the values a > 2/n ~ 1, a > ao the only critical points are those of (|3.2]), (6.15) and ( |6.17 ). The 
corresponding phase portraits are presented in figure 4a. All expanding models start at some finite moment 
of the E- frame comoving time from the point (x = — \Jn(n — 1), w — 1/2). For the closed models they 
terminate at a point (x — — y/n(n — 1), w = —1/2) at another finite time moment. For the op en m odels 
the trajectories approach to the flat spacetime as t — > +oo with the asymptotic behaviour (6.20) and 
(S.21). The spatially-flat models with w — 1/2 (—1/2) are early (late)-time attractors for the open and 



closed models. All models are noninflationary. 



2. When a > 2/n — 1, < a < ao, critical points are those of (6^), (3J5), ( 6.1 5| ) and ( 3.17 ). First two ones 
lie on the segments w — ±1/2 and correspond to the spatially-flat models. Phase-space diagram of this 
case is plotted in figure 4b. As we see for the closed models there are three qualitatively different classes of 
solutions separated by the separatrixes of the saddle points (|6~5|): (i) models with decreasing dilaton (x < 
0), originating at point (—y/n(n — 1), 1/2) in the finite past and approaching point (—y/n(n — 1), —1/2) 
at the finite future, (ii) models with increasing dilaton (x > 0) originating at point (y/n(n — 1), 1/2) in 
the finite past and approaching point (yn(n— 1), 1/2) at the finite future, (iii) models originating at 
point (—^/n(n — 1), 1/2) in the finite past and approaching point (y/ n(n — 1), 1/2) at the finite future. 
In addition, two special solutions exist with semi-infinite lifetime, corresponding to the separatrixes of the 
saddle ( |6.5|) . The first/second one originate at point (— \Jn(n — 1), 1/2) / (a/26, 1/2) in the finite/infinite 
past and approach point (—a/26, —1/2)/ {yf n(n — 1), —1/2) in the infinite/finite future. 

3. The case a = 1/n, a = corresponds to the radiation-dominated models. Phase-space diagram is plotted 
in figure 4c. It is symmetric under reflection x — > —x. Vertical segments x = represent the GR solutions 
for closed (\w\ < 1/2) and open (\w\ > 1/2) models. As we see only expanding GR solution is the attractor 
for generic solution with varying dilaton. For all solutions dilaton field is monotonic. 



4. When a < 2/n— 1, a > ao the critical points are (6.2), (6.S), ( |6.15 ) and ( |6.17 ). The corresponding phase 



portrait is plotted in figure 4d. For p oints (3.8) one has hi(a) > 1 and it represents open models. We 
have assumed that first of conditions ( 6.11 ) takes place and this point is a sink. In the case of uppe r sig n 
i t is t he late-time attractor for all expanding models with e > and is described by expressions ( 4. 35] ), 
( 4.36 ) in the E- and string frames. The qualitative behavior of th e cl osed models is similar to the case 
1 with difference that now exist the segments of trajectories with (h9) corresponding to the inflationary 
evolution in the E-frame. Note that the closed models contain three stages: early-time dilaton-dominated, 
intermediate inflationary and late-time dilaton dominated ones. The open models also contain stages with 



inflationary evolution lying in the region (4.9) 



In the case a < 2/n — 1, a < ao, hi > 1 the new critical points ( |6.5| ) appear on the segments w — ±1/2 
corresponding to the spatially-flat models. Phase-space diagram is plotted in figure 4e. The qualitative 
behavior of the closed mo dels is similar to the case 2, but now an intermediate inflationary stages exists 
in the regions defined by ( 4J3 ) . For t he o pen m odels there are two classes of solutions separated by the 
unstable separatrixes of the saddles (6J5) and ( 6.17 ). In the case of expanding universe the models of 
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both classes (w > 1/2) start near spatially- flat expanding solution in the finite past with increasing and 
decreasing dilaton for the first and second classes respectively. In addition there are two special solutions 



corr esponding to the unstable separatrixes of points (6.5) and ( 6.17 ). The first one originates at point 



(3.5) in the finite past and approaches point (p.q) with upper sign in infinite future. The second solution 



starts at point (6.17) corresponding to the flat spacetime in finite paste and approaches point ( p.8[ ) in 
infinite future. 

When a < 2/n — 1, a < ao, hi < 1 the critical point ( |6~8| ) corresponds to the closed models and is a 
saddle. The phase portrait for this case is presented in figure 4f. There are two classes of expanding 
open models (w > 1/2) separated by unstable separatrix of point ( |6.17 ). The point ( |6.5[ ) is the late-time 



attractor for both these classes. The trajectories of the first/second class start at point (—^n(n — 1), 1/2) 



/ (y/n(n — 1), 1/2) in the finite past and terminate at point ( |6.5| ) in the infinite future. After the initial 
stage of dilaton-dominated expansion these models entry to the inflationary stage. The separatrixes of 
the saddles ( |6.8[ ) split phase space of closed models (\w\ < 1/2) into seven invariant subspaces. The 
trajectories above the stable separatrixes of the point ( |6.8] ) originate at points (±-^/n(n — 1), 1/2) in the 
finite past and approach point ( |6.5| ) at the infinite future. Hence these models have semi-infinite lifetime. 
The same is the case for the closed models below the unstable separatrixes of the saddle (^^) with lower 
sign. The qualitative behaviour of the other closed models is the same as for corresponding ones in the 
cases 2 or 5. They have finite lifetime. 



In the case a = a^\Jl — l/nb (hi = 1) critical points (6.5) and ( pTq ) coincide and correspond to the 
spatially-flat models. As it considered above the critical point is an equilibrium state with one stable node 
sector corresponding to the open models and with two saddle sectors corresponding to the closed models. 
The corresponding phase-space diagram for the open/closed models is qualitatively equivalent to that of 
figure 4 f/c. 

For the nongravitational source with a = 1 one has b = 0. In this case all points of the segments 
w = ±1, < y/n(n — 1) are critical points. The corresponding phase-space diagram on the rectangle 
\x\ < y/n(n — 1), |w| < 1 looks like of that in figure 1. 



In the phase-space diagrams of figure 4 the point (yjn(n — 1), —1/2) corresponds to the pre-big bang gravi- 
dilaton models. As we see this solution is late-time attractor for all contracting models in the cases of 1 and 4 
(figures 4a and 4d) and for the part of models in the other cases of figure 4. 

We have considered tree-level phase portraits for a > 0. The corresponding phase space diagrams for 
a < can be obtained by taking into account the invariance of the dynamical system with respect to the 
transformations x — > — x, a — ► —a. 



7 Qualitative analysis in general case 

Having discussed the cases of spatially-flat and tree-level curved models we are now in a position to consider the 
general case of spatially curved cosmological models without specifying dilaton couplings in the string effective 



action. In the E- frame these models are described by the dynamical system ([h2|). The critical points of this 
system are the points of the phase space (4>,x, h) with coordinates 

(0 O ,O,±1), «(0 O ) = O (7.1) 

where upper/lower sign corresponds to the expanding/contracting models. We shall consider the first of these 
cases. The corresponding results for the contracting models can be obtained by means of transformation ( |4.3| ). 
Note that the points (^J) are GR solutions for k = 0. The constant dilaton solutions for k — 1 are presented 
by segments (j) — 0o, x — 0, \h\ < 1, and for k = — 1 by rays (f> = (f>o, x = 0, \h\ > 1. 

Near the critical point the equation for h decouples from the other two ones and the corresponding 

eigenvalues are equal to 

Ai, 2 = -h ± ■sjkl + n(n - l)a' /2, A 3 =n(l + a)-2 (7.2) 
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The eigenvalues A1.2 define the character of the critical point for the trajectories with h = 1, corresponding to 
the spatially- flat models. The different possible variants of this case have been considered in section 5. As it 
follows from (7.2) for the models k = ±1 the critical point is unstable when a > 2/n — 1. The corresponding GR 
solutions tend to the spatially- flat solution (point ( |7.l| ) with the upper sign) as t — > —00. In the case a < 2/n — 1 
the GR solutions with k = ±1 tend to the solution ( [7.1] ) in the limit t — > +00 (inflationary models). When 
additionally ag < (and therefore ReA^ < 0), the critical point is stable and the solutions of string cosmology 
with variable dilaton tend to the solution (7.1) too. Note that for point (<^0j0, — 1) one has Ai > and this 
point is always unstable. 

Now we turn to the investigation of the behaviour of the trajectories at infinity of the phase space. For 
this it is convenient to compactify the phase space. In the case of the closed models it is sufficient to use the 
mapping (4.15). The corresponding dynamical system differs from ( |4.2| ) by the first of the equations, which now 
has the form 

% = xy{l - y) (7.3) 
and the phase space is the rectangular parallelepiped 



0<y<l, M < y/n(n-l), \h\<l 



(7.4) 



The boundaries of this parallelepiped are invariant subspaces. The trajectories on the boundaries h = 1,-1 
correspond to the spatially-flat expanding and contracting models, respectively, and the trajectories on the 
boundaries x — ±y^n(n — 1) repr esent the pure gravi-dilaton models. For the latter case the equation of the 
phase trajectories have the form ( 4.14 ). 

For \h\ — 1 the critical points on the boundary of fl7.4|) are the points (y,x,h) with coordinates (a = a(y)) 



(j,±\/n(n-l),ho\ > 3 = 0,1, h 
(j,a(j)/2b,h ), \a(j)\<a 



1 (7.5) 

(7.6) 

When ho = 1 (the case ho — — 1 can be obtained by transformation ( [4.3[)) the eigenvalues Aj.,2 for these points 
are the same as in the case of spatially flat models (formulas (5.21), ( |5.22| )), and 



for the point (|7.5|), and 



A 3 = 2(n - 1) 



A 3 = 2n(n - 1)6 [a 2 (j)/a 2 + l/nb-l] 



(7.7) 



(7.8) 



for the points (7.E). Thus, for the models with k = ±1 points (7J3) with ho — 1 are always unstable. Points 
(J, = f( — iy\/n(n — 1), ±1) are saddles always and correspond to an infinite value of the E-frame comoving time. 
Points (j,±(-iy y/An — 1),±1) with upper/lower sign are unstable/stable nodes when (— Vfa.(j) < ao and 
saddles otherwise. Th ey c orrespond to a finite value of the comoving time coordinate and are singular. 
The critical point (|7.6|) with ho = 1/ — 1 is a stable/unstable node when 



ao 



y/i - l/nb < (-l) j a(j) < 0, 



(7.9) 



and a saddle otherwise. Note that the left of these conditions is possible for a < 2/n — 1 only. For the expanding 
models point (7.6) corresponds to a finite/infinite value of the comoving time when (— l) J a(j) > / < 0. By 
taking into account ( [7^ ) we conclude that for k = ±1 all critical points with h = 1 are unstable when a > 2/n— 1. 
In the case a < 2/n — 1 there are critical points laying on the boundaries y = 0, 1 of parallelepiped ( |7.4| ): 

-0,1 



(j, ±y/(n-l)(n-l/b), ±hi(a(j)j) , j = 0, 



(7.10) 



where the function h\{a) is defined by the relation ( |6.8| ). These points correspond to the tree- level ones ( |6.8| ). 
For the points ( 7.10[ ) the eigenvalues A2,3 are determined by the relations ( |6.9| ), and 



Ai = ±(-l)V(n-l)(n-l/&) 



(7.11) 
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When h\ < 1 the points ( 7.10| ) correspond to the closed models and are unstabl e, as the eigenvalues A 2 and 



A3 have different signs. When h\ > 1 the solutions corresponding to the points ( [7.f ) describe open models. 



For the lower/upper sign point (7.f0) with j = is stable/unstable source when a(y = 0) < 0. In the case of 
j = 1 and a(y = 1) > the point ( [7.10 ) is stable/unstable source for the upper/lower sign. For other cases it 
is a saddle. Note that the stable sources correspond to the expanding models with negative spatial curvature. 
In the case of expanding/contracting models the value of the E-frame comoving time for the point (7.1C) is 
t = +00/ — 00 when ±(— iy < j > 0, and this point corresponds to the finite time moment otherwise. For the 
latter case the point ( 7.1 0| ) presents a singularity for cosmological solutions. 



As for the open models \h\ > 1, to compactify phase space along side with (4.15) we shall to map 



z=l/h, \z\<l (7.12) 
By introducing the independent variable T, dT = hdr, the corresponding dynamical system can be written as 



J| = xyz(l-y), ^ = [n(n- 1) - x 2 ] [az/2 - bx] 
^ = z(l-z 2 )[(n-l)(nb-l)-bx 2 ] 



(7.13) 



As in above we shall consider the expanding models (the case of co ntracting models can be obtained by trans- 
formation x — > —x, z — > — z ). The phase space of the system ( 7.1 3| ) is the parallelepiped 



< y < 1, \x\ < y/n{n- 1), < z < 1 



(7.14) 



and the critical points lie on the boundaries z = 1 and z — 0. For the first of these cases (h = 1) the critical 
points and their stability are the same as in the previous case and correspond to the spatially- flat models. Thus 



we shall consider the second case. At z = (h = 00) the solutions of the system (7.13) are the segments z = 0, 
y = constant. All points (y, x, z) of the segments 



(y,x = 0,z = 0), < y < 1 



and 



(y, x = ±\Jn(n - 1), z = 0) 
are critical. The eigenvalues corresponding to the axis z are equal to 

A 3 = (n-l)(n6-l) 



for the points ( |7.15| ), and 



A, 



(7.15) 
(7.16) 

(7.17) 
(7.18) 



for the points (7.16). In the second case the critical points are unstable. The points (7.15) are stable when 
a > 2/n — 1 and unstable otherwise. By taking into account that for the first of these cases the critical points 
for z = 1 are unstable, we conclude that all models with k — — 1 finish evolution on the segment x — 0, z = 0. 
Near the corresponding critical point the trajectories of the dynamical system enter into this point lying in the 
plane (x, z). In this plane the corresponding critical point is an stable node with the eigenvalues 



A 2 = -nb(n - 1), A 3 = (n - l)(nb - 1) 



(7.19) 



In the case a < 2/n — 1 the final point of the evolution for the models with k = — 1 are the points on the 
boundary z = 1, corresponding to the models with k — 0, and the points ( 7.1 C| ) under the conditions specified 
above. 

On the basis of the analysis carried out above it is possible to construct the phase diagrams of the cosmological 
models for arbitrary function a(4>), defined by dilaton coupling functions in the string effective Lagrangian. For 
given a the position and character of the critical points are determined by the zeroes of a(<p) and by the boundary 
values a(±oo). In the parameter space (a, a(y = j)), j = 0, 1 the various qualitatively different regions are 
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separated by curves ( 6.23| ) with a — a(j). As in the tree-lev el ca se, to plot the closed and open models on the 



same phase space diagram we shall use the variable w from ( 4.1 6| ) . In terms of the variables (y, x, w) the phase 



space of the models with nonncgative energy density is the parallelepiped 



0<y<l, |x| < yjn(n-l), \w\<l (7.20) 

The trajectories lying on the sections w — ±1/2 represent invariant subspaces and correspond to the spatially- 
flat models. For expanding solutions the corresponding two-dimensional phase portraits are presented in figures 
2 and 3. The trajectories in the invariant subspaces x — ±^Jn{n — 1) correspond to the pure gravi-dilaton 
solutions presented in figure 1 for increasing dilaton (upper sign). The phase portraits on the boundaries y = 0, 1 
are the same as in the tree-level case with a = a(j) = const (see figure 4). The spatially-flat sections w = ±1/2 
split phase-space into three invariant subspaces corresponding to closed (| w| < 1/2) and expanding (w > 1/2) 
and contracting (w < —1/2) open models. If (ft = (fto is a zero of the function a(<f>) then the dynamical system has 
solutions with constant dilaton. These solutions are presented by the vertical segments {(ft = (j)Q, x = 0, \w\ < 1} 
and correspond to the GR solutions. The horizontal segment x = y/n(n — 1), w = — 1/2 presents spatially-flat 
pure gravi-dilaton contracting solution with increasing dilaton. At tree-level the corresponding string frame 
solution describes a superinflationary type expansion and lies in the basis of pre-big bang inflationary model 
ipHI] , [p4| . As it follows from the above analysis this solution is the late-time attractor for the spatially curved 
models with additional sources when a((ft = +oo) > — a®. 

Let us consider two qualitatively different regions of the parameter a separately. 

7.1 Qualitative evolution: k = ±1, a > 2/n — 1 



In this case nb < 1 and as it follows from ( 4.2) the functions h[r) and w(r) are decreasing/increasing for 
the closed/open models. There are no inflationary segments of trajectories. As possible attractors one has the 
following critical points of the phase space (y, x, w): 

M,±l) (7.21) 
(0,Tv/n("-l)-Tl/2), a(0) < a (7.22) 



(l,±Vn(n-l),=Fl/2), a(l) > -oo (7.23) 
where upper/lower sign corresponds to the late/early-time attractors. The other critical points are saddles and 



hence are unstable. The contracting/expanding open models start/finish their evolution at the points (7.21) 
with lower/upper sign, where y is determined by the initial conditions. The initial stage for the expanding 
models is dilaton dominated and corresponding trajectories originate close to the pu re gr avi-dilaton spatially- 



flat solutions. During the further expansion the open models approach the point ( 7.21 ) with u p per s ign at 
the infinite future. In a neighborhood of this point the solution have the asym ptotic form (|6.20|) , ( 6.21 ) with 



a = a(y) and (ft tends to a constant. Assuming the existence of the zeroes (7T) with a' ((fto) < for the closed 
models the field (ft will initially evolve towards the value (ft = (fto. During the further evolution h becomes zero 
at some finite time moment and the universe enters into the contracting phase. In this stage the friction term 
in the scalar field equation changes sign taking away the scalar field from the critical point. Hence in this case 
the GR solutions are not late/early-time attractors for the general solution with varying dilaton. 

Let us consider the qualitative evolution of the spatially curved models on the phase-space (y, x, w) for 
different boundary values a(j), j = 0, 1. 



1. a(j) > eta, j = 0, 1. The critical points are possible zeroes of the function a(cj>), ([74|), the points ( 7.1 5| ) 



and ( 7.1 6| ) . All expanding models with k = ±1 and e > 0, apart from the possible trajectories with 



constant dilaton, originate at point (7.23) with lower sign in the finite past of the E-frame comoving time 



and approach point (7.23) with upper sign in the finite future for the closed models, and point (7.21) with 



upper sign in the infinite future for the open models. In the latter case the limitary values of dilaton are 
determined by the initial conditions. The pre-big bang gravi-dilaton solution x — ^Jn(n — 1), w = —1/2 
is the attractor for all contracting models with k = ±1. The corresponding 3D phase-space diagram is 
plotted in figure 5a, assuming that the function a((ft) has no zeroes. In this diagram the phase portraits 
on the boundaries y = 0, 1 are similar to that of figure 4a. 
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a(0) > Ooj < a(l) < ot$. In addition to the critical points of the previous case, the new unstable critical 
points ( |7.6| ) with j = 1 appear for spatially-flat models. The qualitative behaviour of the trajectories with 
k = ±1 (w ±1/2) is the same as in the case 1. When the function a(4>) has no zero the phase portrait 
on the section w = 1/2 is similar to that of figure 2e. Phase portrait on the boundary y = 1 is equivalent 
to that of figure 4b. 



3. a(0) > ao, — ao < a(l) < 0. Now the critical point (7.6) is unstable for the spatially-flat models. The 



qualitative picture for the spatially-curved models is the same as in the case 1. 

a(0) > ao, Q;(l) < — ao- For the spatially-flat models there are no stable critical points on the phase space 
boundary. The boundary rectangles for the spatially-flat models are limitary cycles for the spatially- 
curved solutions. All expanding models with k — ±1 originate close to the rectangle with w = 1/2 in the 
infinite past and approach to the rectangle with w = —1/2 in the infinite future for the closed models. The 
late/early-time behaviour of the expanding/contracting open models is the same as in the previous cases. 
The corresponding phase diagram is plotted in figure 5b. The phase portrait in the invariant subspace 
w = 1/2 is similar to that of figure 3a. 



< a(0) < ao, a(l) > ao- The critical point (7.6) with j = corresponds to the spatially-flat models and 
is a saddle. The expanding open models originate at points (7.22) and (7.23) with lower signs in the finite 
past and approach point (7.21) at the infinite future. The corresponding limitary value of y depends on 
the initial conditions. In the simple case when a((f>) > for the closed models there are three classes of 



solutions: (i/ii) models originating at point (7.23) with lower sign in the weak coupling region at the finite 
past and approaching point ( [7.22 )/( 7.23] ) with upper sign at the finite future, (iii) models originating at 
point (7.22) with lower sign in the finite past and approaching point ( [7.23 ) with upper sign. Note that 
(dx/ dr) x= Q — n(n — l)a((f)) /2 > and therefore the trajectories originating from (7.22) with lower sign 
can not terminate at point ( [7.22 ) with upper sign. The early and late-time asymptotics of these solutions 
are pure gravi-dilaton spatially-flat models. For the trajectories of the class (ii) both these asymptotics 
lie in the strong coupling region. 



< a(j) < a , j = 0,1. In addition to the points of the previous case the new critical point, (7.22), 
j = 1 appears. But this point corresponds to the spatially-flat models and the qualitative features of the 
spatially curved models remain the same as in the previous case. 

< a(0) < ao, — ao < a(l) < 0. The function a (0) has a zero with negative derivative. The corresponding 
critical points are unstable. The qualitative behavior of the trajectories for the spatially curved models 
with variable dilaton remains the same as in the case 5. For the universe with positive curvature in 



addition to the models of the case 5 there is also a class of solutions originating at point (7.22) with lower 
sign in the finite past and approaching point (7.22) with upper sign in the finite figure.. The corresponding 
phase diagram is plotted in figure 6a. 

< a(0) < ao, a(l) < —ao. For the models with k — ±1 all critical points are saddles, except the 
points (7.22) and (7.21). The first of these points with lower/upper sign is early/late-time attractor for 



expanding/contracting models. All closed models originate in the weak coupling region in the finite past 
and terminate in the same region at the finite future. 

Radiation-dominated models: a = l/n. The corresponding phase space diagram is presented in figure 6b. 
The vertical segments {y = const, x — 0, w) are solutions of the dynamical system and correspond to the 
models with constant dilaton. They split the phase-space into two invariant subspaces with increasing (x > 
0) and decreasing (x < 0) dilaton. The solutions of the first/second class originate at point ( [7.22 ) / (7.23) 
with lower sign in the finite past and approach point (7.21) in the infinite future for the open models and 
(f?^/(|7\2|) with upper sign in the finite future for the closed models. 



7.2 Qualitative evolution: k = ±1, a < 2/n — 1 



Now there are segments of the trajectories with accelerated evolution lying in the region (4^). When a < 
the point (7T) with upper/lower sign is late/early-time attractor for the models with varying dilaton. As in 
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the case of the spatially-flat models this provides an efficient mechanism for dilaton stabilization during the 
expansion of the universe and is the basis of the Damour-Polyakov mechanism. The other critical points lie on 
the boundary of the phase-space (y,x,w). In addition to the points ( 7.22j ) and ( [7.23 ) as possible attractors at 
the phase-space infinity one has the following critical points 



(j, ±a(j)/2b, ±1/2) 



-a 



y/l-l/nb < (-l)Mj) < 



j, t(-1)V("-1)(»- 1/6). 



T(-i)Mj) 



aoy/l — 1/nb 



(-l) j a(j) < -a oy /l - 1/nb 



(7.24) 



(7.25) 



where upper/lower signs correspond to the late/early-time attractors. Other critical points (including (7.21)) 



are saddles. The points (7.25) correspond to the open models. In a neighborhood of these points the solution 
has the form ( |6.20| ), fl6.21|) with a = a(j). When j — dilaton is in the weak c oupli ng r egion and the 
corresponding asymptotics for string frame scale factor and dilaton are in the form ( 4.36 ) or ( 4.38| ). Let us 
consider qualitatively different phase diagrams. 



a(j) > do, j = 0,1. On the boundary of the parallelepiped ( [7.20 ) one has late/early-time attractors 
( [7. 22 ) and (7.23). The corresponding phase diagram is plotted in figure 7a assuming that a((f>) > 0. In 
this case there are no critical points at the finite part of the phase-space. The expanding models originate 
at point (7.23) with lower sign in the finite past and approach (i) point (7.23) with upper sign at the 
finite future for k = 1 and (ii) point ( 7.25| ) with upper sign at the infinite future for k = — 1. In the 
case of models with positive spatial curvature h becomes zero at some finite time moment and then the 
universe enters into the contracting phase. Note that when a < 2/n—l the function h(r) is not monotonic 
and in general case the expansion-contraction transition can take place several times. After the initial 
dilaton-dominated decelerating expansion the models enter into the inflationary stage which corresponds 
to dh/dt > 0. For the closed models after some amount of inflation the trajectories enter again into the 
deflationary phase. After inflation the open models enter into the oscillatory regime in a neighborhood of 
point ( 7.25| ) with upper sign. In addition to these models there are also special solutions corresponding 
to the unstable separatrixes of the points ( 7.21 ) with upper sign (see figure 7a). In the case a(0) > ao, 
yl — 1/nb < a(l)/aa < 1 the qualitative features of the phase diagram for the spatially curved models 
is the same as in the figure 7a (note that this is not the case for the spatially-flat models). 



2. When a(0) > ao, < a(l) < ao\/l — 1/nb the critical point (7.25) is a saddle and corresponds to the 
closed models. Now at the infinity of the phase-space one has attractors (7.23) and (7.24) with j = 1. 
The qualitative behaviour of the open models is the same as in the previous case, with difference that 
all expanding models approach point (7.24) with j = 1. For the models with k = 1 in addition to the 



trajectories considered above there are models approaching the point (7.24) with j = 1 at the infinite 
future. These models approach to the spatially-flat ones with the infla tiona ry type final stage. There are 
also special solutions corresponding to the stable separatrix of saddle ( 7.25 ), j = 1. 



a(0) > ao, —yl — 1/nb < a(l)/ao < 0. The function a(<fi) has a zero with a'(<f>o) < and the point 
( [7.l| ) with upper sign is an attractor. Phase diagram for this case is plotted in figure 7b. The segment 
x = 0, y = y((f>o) presents the corresponding GR solutions. The expanding models with k = —1 originate 
at point ( 7.23| ) with lower sign in the finite past and approach point (7.1) with upper sign at infinite 
future. There are also special solutions corresponding to the unstable separatrixes of the saddles (7.21). 
For these models one has efficient dilaton stabilization. In the case of the universe with positive spatial 
curvature there are four classes of solutions having as early and late-time attractors the points (i) ( 7.23| ) 
with lower sign and (7.23) with upper sign, (ii) (7.23) with lower sign and ( [7.1[ ) with upper sign, (iii) point 
(7.1) with lower sign and (7.1) with upper sign, (iv) point (7.1) with lower sign and ( [7.23 ) with upper 
sign. Dilaton stabilization by means of Damour-Polyakov mechanism takes place for the models (iii) only. 
In the case of the models (i) the initial and final stages of the cosmological evolution correspond to the 
dilaton dominated decelerating expansion and contraction respectively, whereas the intermediate stage is 
perfect fluid dominated with inflationary type evolution. The models of (ii) and (iv) have semi-infinite 
lifetime. In the case (iii) the universe originate in the infinite past at point (7T) with lower sign. After 
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initial contraction the universe enters into the expanding stage and approach point (7.1) with upper sign 
as a result of inflation. For parameters region a(0) > do, —1 < a(l)/ao < — y/1 — 1/nb the qualitative 
features of the phase-space diagram for the spatially curved models is similar to that of figure 7b. 

a(0) > olq, cs.[l) < — otQ. There are no attractors on the phase-space boundary. The boundary segments 
of the rectangles w — ±1/2 are early (for upper sign) and late-time (for lower sign) limitary cycles. The 
qualit ative bahaviour of the models is the same as in the previous case with difference that instead of nodes 
( 7.23 ) now one has limitary cycles. The corresponding phase-diagram is plotted in figure 8a, assuming a 
single zero for the function a((j>). 



5. yjl - 1/nb < a(0)/a Q < 1, a(l) > Qj Q. In the case of positive function a(<p) for the expanding open models 
the only late-time attractor is point (7.25) with j = 1. There are two classes of such models originating at 
points ( [7.22 ) and ( 7.23| ) with lower signs in the finite past. In addition special solutions exist corresponding 

with j = 0. The closed models 



to the unstable separatrixes of points ( 7.21 ) with uppe r sig n and (6. 

can be divided into three classes: originating at point (7.23) with lower sign and approaching points (i) 



(7.22) and (ii) ( 7.23|) wi th upper signs, and trajectories originating at point (7.22) with lower sign and 
approaching point (7.23) with upper sign. As in the case 5 for a > 2/n — 1 in the case of the positive 
function a(<p) the trajectories originating at point (7.22) with lower sign can not terminate at point ( 7.22| ) 
with upper sign. 



^1 - 1/nb < a(0)/a o < 1, < a(l) < a y/l - 1/nb. As an attractors one has Q7.22D , ( [7.23[ ) and fl7.24| ) 
with j = 1. The qualitative features for the open universe is similar to those of the previous case with 



only difference that now instead of point (7.25), j — 1 one has point (7.24), j = 1 as late-time attractor 
for the expanding models. In the case of the closed universe two new classes of solutions appears (and 
the solutions obtained from these ones by time inversion) in addition to previous case. They originate 
at points ( 7.22j ) and ( 7.23| ) with lower signs in the finite past and approach the point (7.24), j = 1 at 
the infinite future. The corresponding phase-space diagram is depicted in figure 8b assuming a positive 
function a (</>). 



i/l — 1/nb < a(0)/ao < 1, — aoyl — 1/nb < a(l) < 0. The function a((j>) has a zero with negative 
derivative. The phase-space diagram for the case of single zero is presented in figure 9a. Critical points 



( |7.l|) are e arly (lower sign) and late-time (upper sign) attractors. The other attractors are points (|7.22| ) 
and (7.22). The expanding open models originate at these points with lower sign in the finite past and 
approach point (7.1) with upper sign at the infinite future. The evolution consists two stages: initial 
decelerated dilaton-dominated and final inflationary type one. As a result of the cosmological expansion 
the dilaton variations are strongly suppressed. There are also special open solutions corresponding to 
t he u nstabl e se paratrixes of points ( 7.21 ) with upper sign. The closed models originate at points ( 7. 21 ) , 
(7.23) and (7.1) with lower signs and approach these points with upper signs. As a consequence one has 
nine qualitatively different types of the closed models. There are also special solutions corresponding 
to the separatrixes of the saddle point. For the models having the point (7.1) with upper sign as a 
late-time attractor the dilaton is efficiently stabilized. Final stage of these models present inflationary 
type expansion dominated by barotropic perfect fluid with a < 2/n — 1. There are also models with 
intermediate inflation. As we see from figure 9a, depending on the initial conditions the both types of 
the inflationary evolution: standard and pre-big bang, can be realized within the framework of the same 
model and correspond to the trajectories having as a late-time attractor points ( |7.l| ) and (7.22) with upper 
signs. 



— ao\/l — 1/nb < a(0) < 0, a(l) > olq. The function a((f>) has a zero with a'(</>o) > 0. The corresponding 
critical point is a saddle and the solutions with constant dilaton are unstable. In figure 9b we have plotted 
phase space diagram assuming a single zero for the function a ((f)). All attractors lie on the phase-space 
boundary. These are points ( |7.22D , ( [7.23j ), ( |7.24| ) with j = and ( |7.25| ) with j = 1. The latter point 
is a late-time attractor for the expanding open models. These models can be divided into three classes 



with qualitatively different behaviour originating at points (7.22), (7.23) with lower sig ns an d ( 7.24|) with 
j = 0. All closed models except the solution with constant dilaton originate at points ( 7.22 ), ( 7.23| ) with 
lower signs in the finite past and approach the corresponding points with upper signs at the finite future. 
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In the case — ao\A — < < 0' \A — l/nb < g(l)/ ag < 1 the phase-space diag ram is similar to 
that of previous case with difference that now one has point ( 7.24 ) with j = 1 instead of ( 7.25 ) with j = 1. 
The latter corresponds to the closed models and is a saddle. For the closed universe in addition to the 



models of previous case new class of solutions appears having as a late-time attractor the point ( 7.24 ) 
j = 1. This point correspond to the spatially-flat models. 



The phase diagrams for the another combinations of the values of a(0) and a(l) are qualitatively equivalent 
to the ones considered above or can be obtained from those by taking into account the invariance of the 
dynamical system under the transformation (4.4). 



8 Conclusion 

We considered the evolution of a homogeneous and isotropic universe within the framework of the string effective 
gravity with the string-loop modifications of the dilaton couplings, (note that we have investigated the most 
general case of coupling functions and our analysis is valid for any _D-dimensional scalar-tensor theory. For the 
barotropic perfect fluid as a nongravitational source th e co rresponding set of equat ions can be presented in the 



form of the third-order autonomous dynamical system (h2) by using the variables (4J). The model is sp ecifie d 



by the barotropic index a and by the function a{4>). The latter is determined by dilaton couplings (see ( 3.1 2| ) ) 



and is a constant at tree-level. As an example of an additional source we consider form-fields, which naturally 
arise in supergravities deriving from string theories. 

For the pure gravi-dilaton, radiation-dominated and stiff fluid models the solutions to the dynamical system 



(4.2) can be found in terms of integrals containing combinations of the dilaton coupling functions. In the case 



of constant a (w hich include tree-level models as a particular case ) the re are two power-law special solutions, 



( 4.24 ) and ( 4.34 ) with the E- frame time dependences ( 4.25 ) and ( 4.35 ), respectively. In addition to the pure 



gravi-dilaton they are early or late-time attractors for the models with varying a. The first of these solutions 
corresponds to the spatially-flat models and describes an extended inflation (in the E-frame) in the case of 



(4.26). In this case the second special solution describes the universe with positive spatial curvature. 

For the general case of a and function «(</>) to obtain the generic features of the cosmological evolution we 
have used the dynamical systems methods. We start with spatially-flat models in section 5 when dynamical 
system is reduced to the second-order one. The only critical points at the finite part of the phase space are those 
corresponding to the zeroes of «(</>) (for radiation a — and all points of the segment x — are critical). These 
points are saddles when a'(</>o) = a' Q > 0, centers for a' Q < 0, a = 1, and late-time attractors (for expanding 
universe) in otherwise. In the last case we have an efficient dilaton stabilization as a result of the cosmological 
expansion (Damour-Polyakov mechanism). Note that in the case of the spatially- flat models this mechanism 
works for any a < 1, despite the fact that the stabilization efficiency depends on a. The most favorable case is 
a = — 1. The stability of the critical points at the phase space infinity is investigated and various qualitatively 
phase portraits are presented in figures 2 and 3. Another case when the dynamical system is reduced to the 
second-order one are models with constant a. The phase space analysis of these models is given in section 6 
and results are presented in figure 4. Now dilaton can not be stabilized at finite values. 

The general case of the spatially curved models with arbitrary function a{4>) is considered in section 7. 
The corresponding phase space diagrams are three-dimensional. Their qualitative structure are specified by the 
zeroes of a(4>), limitary values a{4> = ±oo) and barotropic index a. The qualitatively different regions of the 
latter are separated by a = 2/n — 1. The critical points corresponding to the zeroes of a(</>) with a' < 0, a < 1 
are attractors for the models with curved space when a < 2/n — I only. Another attractors tie at the infinity of 
the phase space and correspond to the gravi-dilaton or special solutions (4.24) and (4.34). For these values of a 



there arc inflationary segments of trajectories lying in the region (4.9). For the models having the zeroes of a((f>) 
as a late-time attractor the final stage of the evolution is an inflationary type expansion with efficient dilaton 
stabilization. As dilaton variations in this stage are small it corresponds to the inflation in the string frame 
as well. There is also a class of trajectories with pre-big bang gravi-dilaton solution as a late-time asymptotic. 
Note that depending on the initial conditions the both types of the inflationary evolution: standard and pre-big 
bang, can be realized within the framework of the same model. For various values of parameters the different 
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scenarios of the cosmological evolution are described in subsection 7.1 for a > 2/n — 1 and 7.2 for a < 2/n — 1. 
The corresponding phase space diagrams are presented in figures 5-9. 

In this paper we have considered the evolution of the cosmological models in the E-frame. They have 
singularities at a finite E-frame comoving time. In this higher curvature regime higher derivative terms due to 
string corrections become important. The string frame solutions can be obtained by transformations 

t s = J n E (<p)dt, Re = n E (<p)R (8.1) 

with frame comoving time t s . By choosing dilaton coupling function one could construct the metric in the 
string frame which is singularity free over an infinite interval of the string comoving time. Analogous models 
for the Jordan frame in the scalar-tensor theories were considered in J83| , |S4| , where constraints are specified 
which lead to a singularity free evolution in the Jordan frame. However as it have been shown in fl67f the 
conformal transformations alone cannot completely remove the singularities from cosmological solutions in the 
scalar-tensor models. In particular, gravitons follow geodesies in the E-frame, which are incomplete and have a 
singularity at a finite E-frame comoving time. Hence gravitational waves see the singularity. 
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Figure captions 

Figure 1. Phase portrait for the gravi-dilaton models. The phase space is mapped onto rectangle (y,w), 
< y < 1, < 1 according to ( 4.15| ) and ( 4.16| ). The horizontal segments w = ±1/2 correspond to the 
spatially-flat models. The trajectories lying in region \w\ > / < 1/2 describe the universe with negative/positive 
spatial curvature. 



Figure 2. Phase portraits for the spatially-flat models in the phase space (y, x), < y < 1, |x| < ^/n(n — 1. 
(a) tree-level models with < a < ao, (b) tree-level models with — ao < a < 0, (c) radiation-dominated models 
a = 0, (d) a((f>) > 0, < a(±oo) < ao, (e) ce(4>) > 0, a(— oo) > ao, < a(+oo) < ao, (f) — a o < a(— oo) < 0, 
< a(+oo) < ao with a single zero. 

Figure 3. The same as in figure 2. (a) a(— oo) > a , a(+oo) < — a with a single zero, (b) < ce(— oo) < a , 
a(+oo) < — ao with a single zero as an improper node, (c) < ce(— oo) < ce Q , —ao < a(+oo) < 0, with a single 
zero as a node, (d) — a < a(±oo) < with two zeroes, (e) stiff fluid (a = 1), a(±oo) < with two zeroes, (f) 
stiff fluid, a{— oo) > ao, a(+oo) = 0. 

Figure 4. Phase portraits for the models with constant a compactified onto rectangle (x,w), \x\ < 



\Jn(n — 1, \w\ < 1. Vertical segments x — ±-y/n(n — 1) correspond to the gravi-dilaton models. The tra- 
jectories lying in region \w\ > / < 1/2 describe the universe with negative/positive spatial curvature, (a) 
a > 2/n — 1, a > ao, (b) a > 2/n — 1, < a < a , (c) radiation-dominated models, a = 1/n, a = 0, (d) 
a < 2/n — 1, a>ao, (e) a < 2/n — 1, a<ao, hi > 1, (f) a < 2/n — 1, a <ao, hi < 1. 

Figure 5. Phase space diagrams for the spatially curved models with a > 2/n — 1 compactified onto 
parallelepiped (y,x,w), < y < 1, \x\ < y / n(n — 1, \w\ < 1. Sections w = ±1/2 correspond to the spatially-flat 
models. The trajectories lying in region \w\ > / < 1/2 describe the universe with negative/positive spatial 
curvature, (a) a(y = 0, 1) > ao, a(y) > 0, (b) a(0) > ao, a(l) < — ao. 

Figure 6. The same as in figure 6. (a) < a(0) < ao, — ao < a(l) < 0, (b) radiation-dominated models, 
a = 1/n, a = 0. 

Figure 7. The same as in figure 5 for models with a < 2/n — 1. (a) a(y = 0,1) > ao, a(y) > 0, (b) 
a(0) > ao, - l/nb < a(l)/a < 0. 



Figure 8. The same as in figure 7. (a) a(0) > ao, a(l) < — ao, with a single zero for a(y), (b) — l/nb < 
a(0)/ao < 1, a(l) > ao, with positive function a(y). 



Figure 9. The same as in figure 7. (a) yl — l/nb < a(0)/ao < 1, — yl — l/nb < a(l)/ao < 0, with a 
single zero for a{y), (b) — y 7 ! — l/nb < a(0)/a o < 0, a(l) > a , with a single zero for a(y). 
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